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The  purpose  of  this  study  was  to  determine  if  the 
OSRC  approximation  technique  could  be  applied  to  a  coated 
scatterer  and  achieve  reasonable  accuracy.  My  original 
intentions  were  to  do  much  more  work  with  non-circular 
geometries  than  with  circular  geometries.  However,  as  is 
the  case  so  often  in  research  projects,  unexpected  results 
and  time  limitations  forced  me  to  limit  the  class  of 
problems  considered  to  just  circular  geometries.  The 
results  presented  here  should  be  indicative  of  what  may 
happen  for  a  non-circular  scatterer. 
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Abstract 

The  purpose  of  this  study  was  to  use  On-Surface 
Radiation  Conditions  (OSRC)  to  approximate  the  RCS  of  a  two 
dimensional  dielectric  coated  circular  cylinder  and  compare 
the  results  to  an  exact  solution.  The  comparisons  were  made 
to  determine  a  lower  and/or  upper  bound  on  the  frequency, 
radius  of  curvature  of  the  cylinder  or  any  other  pertinent 
parameters  for  a  valid  approximation. 

The  problem  was • formulated  for  a  circular  cylinder 
of  radius  a  coated  by  a  lossy  dielectric  of  radius  b  and  a 
unit  magnitude  plane  wave  incident  from  the  ^  =  180  degrees 
direction  with  wave  number  k.  The  OSRC  solution  and  the 
eigenfunction  solution  were  nearly  identical;  however,  the 
overall  accuracy  of  the  OSRC  technique  was  shown  to  be 
governed  by  a  complicated  expression  which  was  dependent  on 
the  dielectric  material  parameters  as  well  as  the  kb 
product.  The  solutions  were  compared  for  three  different 
circumstances:  fixed  frequency  and  varying  dielectric  layer 
thickness;  fixed  geometry  and  varying  frequency;  fixed 
frequency  and  geometry  and  bistatic  RCS.  The  OSRC  solution 
was  said  to  be  good  if  the  error  was  less  than  2  db. 

The  OSRC  solution  was  found  to  agree  with  the 
eigenfunction  solution  for  wide  ranges  of  parameters.  In 


general  these  parameters  were  the  thickness  of  the 
dielectric  layer  less  than  0.4  wavelengths  for  TM 
polarization,  the  magnitvfde  of  the  imaginary  part  of  the 
relative  permittivity  greater  than  1.0  for  incident  TM 
polarization  or  1.5  for  incident  TE  polarization,  or  the  kb 
product  "very"  large.  Data  was  presented  which  precisely 
outlined  which  combinations  of  geometry,  dielectric  material 
loss,  polarization  and  frequency  produced  a  valid  OSRC 
approximation . 

These  results  are  consistent  with  the  physical  behavior 
of  the  field  in  the  dielectric  layer.  The  OSRC  technique 
was  shown  to  be  inaccurate  for  higher  order  modes.  Whenever 
a  situation  was  examined  which  would  be  likely  to  excite 
higher  order  modes,  the-  error  in  the  OSRC  solution  was 
large.  For  instance,  TE  polarization  is  much  more  likely 
to  excite  traveling  waves  along  the  surface  of  the 
conducting  cylinder.  These  traveling  waves  are  higher  order 
modes  of  the  field  in  the  dielectric  layer.  Data  was 
presented  which  showed  that  the  error  in  the  OSRC  solution 
was  greater  for  TE  polarization  than  for  TM  polarization, 
for  identical  geometries  and  frequencies,  due  to  the 
contribution  of  these  higher  order  modes. 


THE  RADAR  CROSS  SECTION  OF  A  TWO  DIMENSIONAL 


DIELECTRIC  COATED  CYLINDER  USING  ON-SURFACE 
RADIATION  CONDITIONS 

I.  Introduction 


Background 

Exact  analytical  solutions  for  electromagnetic  scattering 
and  diffraction  phenomena-  are  limited  to  a  small  class  of 
simple  shapes  whose  surfaces  coincide  with  those  of  an 
orthogonal  curvilinear  coordinate  system.  Use  of  a  digital 
computer  makes  it  possible  to  obtain  rigorous  numerical 
solutions  for  electrically  small  structures  of  arbitrary 
shape  using  an  integral  equation  formulation.  However,  due 
to  the  finite  capacity  of  the  computer  this  method  becomes 
untenable  as  the  frequency  increases  (James,  1974:1). 

These  constraints  on  solving  scattering  problems  have  given 
rise  to  many  different  approximation  techniques.  Ruck  and 
others  describe  several  different  approximation  techniques 
for  computing  the  scattered  fields  and  thus  the  radar  cross 
section  (RCS)  of  structures  which  are  not  amenable  to  exact 
solutions.  Most  notable  are  physical  optics  (PC) , 
geometrical  optics  (GO)  and  the  geometrical  theory  of 
diffraction  (GTD)  (Ruck  and  others,  1970:39-105  and  Jost, 
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1988) .  Other  techniques  which  have  been  used  to  improve  on 

.  . 

these  high  frecjuency  approximations  are  physical  theory  of 

diffraction  (PTD)  and  uniform  theory  of  diffraction  (UTD) 
(James,  1974:4  and  Leader,  1982:19).  Many  new  techniques 
and  variations  of  these  techniques  appear  in  the  literature 
every  year. 

A  new  high  frequency  technique,  which  has  been  named  On- 
Surface  Radiation  Conditions  (OSRC) ,  has  been  shown  to 
produce  an  analytical  approximation  for  the  RCS  of  several 
two  dimensional  perfectly  conducting  bodies  (Kriegsmann  and 
others,  1987:153-155).  The  purpose  of  this  study  was  to  use 
OSRCs  to  approximate  the  RCS  of  a  two  dimensional  dielectric 
coated  circular  cylinder  and  compare  the  results  to  an  exact 
solution  to  determine  the  validity  of  the  method. 

Problem  Statement 

The  RCS  of  a  two  dimensional,  perfectly  conducting, 
circular  cylinder  coated  with  a  lossy  dielectric  material 
will  be  found  using  On-Surface  Radiation  Conditions.  This 
approximate  solution  will  be  compared  to  the  exact  solution, 
for  circular  geometries,  to  determine  a  lower  and/or  upper 
bound  on  the  frequency,  radius  of  curvature  of  the  cylinder 
or  any  other  pertinent  parameter  for  a  valid  approximation. 
Any  improvement  in  computational  efficiency  will  also  be 


measured. 


Suitimarv  of  OSRC  Method 


The  OSRC  technique  uses  the  asymptotic  behavior  of  the 
appropriate  Green's  function  to  establish  a  series  of 
differential  operators  that  annihilate  the  first  N  terms  of 
the  asymptotic  series  of  the  scattered  field  (Kriegsmann  and 
others,  1987:153).  '  ' 

Kriegsmann  and  Morawetz  first  derived  the  differential 
operators  to  be  used  in  this  investigation  (Kriegsmann  and 
Morawetz,  1980:371-385);  however,  the  operators  were  applied 
at  some  distance  off  the  surface  of  the  scatterer.  The 
interested  reader  is  referred  to  this  paper  for  the  details 
of  the  derivation.  Kriegsmann  and  others  were  the  first  to 
apply  these  differential  operators  directly  on  the  surface 
of  the  scattering  body  (Kriegsmann  and  others,  1987:153- 
155) .  Their  problem  was  formulated  using  the  free  space 
Green's  function  to  write  a  surface  integral  equation  for 
the  scattered  field.  As  -usual,  both  the  scattered  field  and 
its  first  derivative  with  respect  to  the  surface  normal 
appeared  in  the  integrand.  However,  since  the  scatterer  was 
perfectly  conducting,  the  scattered  field  was  replaced  by 
the  negative  of  the  incident  field  evaluated  on  the  surface 
of  the  scattering  body.  The  problem  was  thus  reduced  to 
finding  an  expression  for  the  normal  derivative  of  the 
scattered  field  on  the  siirface  of  the  scatterer. 

Application  of  a  differential  operator  to  the  asymptotic 
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series  for  the  scattered  field  produces  an  approximation  for 
this  derivative  which  is  accurate  to  order  ^  where 

N  is  the  order  of  the  operator  and  R  is  the  vector  magnitude 
of  the  distance  from  the  surface  of  the  scatterer  to  the 
origin.  The  order  of  the  operator  implies  that  it 
annihilates  the  first  N  terms  of  the  asymptotic  series  for 
the  scattered  field.  The  expression  for  the  normal 
derivative  was  substituted  into  the  integral  equation  and  a 
far  field  expansion  for  the  Green's  function  produced  an 
asymptotic  series  for  the  scattered  field.  The  surface 
electric  current  density  was  found  to  be  proportional  to  the 
sum  of  the  normial  derivatives  of  the  scattered  and  incident 
fields.  The  scattering  cross  section  was  found  to  be 
proportional  to  the  magnitude  squared  of  the  first 
coefficient  of  the  asymptotic  series  representing  the 
scatter  d  field  (Kriegsmann  and  others,  1987:154-155). 

Kriegsmann  and  others  used  the  first  (B^)  and  second  (B^) 
order  operators  to  solve  for  the  scattering  cross  section 
and  surface  currents  of  several  perfectly  conducting 
geometries.  The  first  three  examples  were  for  TM 
polarization  and  the  last  example  was  for  TE  polarization. 

The  bistatic  RCS  for  a  circular  cylinder  was  compared  to 
the  exact  solution  and  found  to  agree  to  within  0.5  dB  for 
and  to  significantly  disagree  for  B^.  The  OSRC  solution, 
using  B^ ,  required  the  evaluation  of  only  two  Bessel 
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functions  regardless  of  the  electrical  size  of  the  cylinder. 
Additionally,  both  OSRC'  solutions  produced  a  non-zero 
surface  current  in  the  shadow  region  of  the  cylinder.  This 
was  shown  to  be  significantly  closer  to  the  exact  solution 
than  the  physical  optics  solution  for  the  surface  currents. 
Next,  the  OSRC  technique,  using  only,  was  applied  to  a 
conducting  strip  to  find  the  monostatic  cross  section.  This 
solution  was  compared  to  a  method  of  moments  (MM)  solution 
and  found  to  be  within  1  dB  for  angles  of  incidence  within 
thirty  degrees  of  broadside,  except  at  nulls.  Disagreement 
at  grazing  angles  was  attributed  to  edge  currents.  Then  the 
OSRC  approximation  was  found  to  agree  within  0.5  dB  with  a 
moment  method  solution  for  the  cross  section  of  a  square 
conducting  cylinder  for  angles  of  incidence  within  45 
degrees  of  broadside.  Finally,  the  surface  current  density 
for  a  circular  cylinder  with  incident  TE  polarization  using 
B^  and  was  compared  to  a  moment  method  solution.  Again 
the  solution  was  found  to  agree  very  well  with  the  MM 
solution  and  produce  a  non-zero  current  in  the  shadow  region 
of  the  cylinder.  In  all  the  cases  presented  above,  an 
analytical  expression  for  the  surface  current  or  scattering 
cross  section  was  derived  and  the  expressions  for  the 
surface  currents  remained  valid  across  the  shadow  boundary 
(Kriegsmann  and  others,-  1-987:155-159). 


Two  approximations  for  the  RCS  using  OSRC  will  be 
presented.  Both  solutions  use  the  two  dimensional  Green's 
function  to  write  an  integral  equation  for  the  scattered 
field.  The  first  method  uses  the  boundary  conditions  at  the 
interface  between  the  dielectric  layer  and  free  space  to 
express  the  integration  in  terms  of  the  incident  field  and 
the  field  in  the  dielectric  material.  Since  the  incident 
field  is  known,  the  OSRC  -technique  will  be  used  to 
approximate  the  field  in  the  dielectric  material.  The 
second  method  uses  the  OSRC  technique  combined  with  the 
boundary  conditions  to  express  the  integral  in  terms  of  the 
incident  field  and  the  field  in  the  dielectric  material. 

Both  methods  use  the  asymptotic  behavior  of  the  Green's 
function  to  derive  a  far  field  expression  for  the  scattered 
field  from  the  integral  equation. 

The  bistatic  RCS  of  a  dielectric  coated  circular  cylinder 
will  be  found  using  the  approach  described  above.  Both  OSRC 
solutions  will  be  compared  to  the  exact  series  solution  for 
this  geometry  (Tang,  1957-; 639-630)  .  Agreement  between  the 
two  solutions  will  be  considered  a  validation  of  the  OSRC 
solution  for  the  ranges  of  parameters  compared. 


Assumptions 

The  discussion  presented  above  is  based  on  several 
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assumptions.  First,  all  fields  are  time  harmonic;  an  exp(- 
jwt)  time  dependence  will  be  assumed  throughout.  Next,  the 
inner  cylinder  is  perfectly  conducting  and  the  dielectric 
material  is  homogeneous.  Finally,  the  surface  charge  and 
current  density  of  the  coating  material  are  assumed  to  be 
zero.  These  assumptions  are  standard  for  this  class  of 
problem. 


Notation 

A  consistent  set  of  notation  was  used  throughout  the 
study.  An  underlined  quantity  is  generally  a  vector 
quantity.  A  primed  (')  coordinate  is  the  integration 
coordinate  and  an  unprimed  coordinate  is  the  observation 
coordinate.  The  cylindrical  coordinate  system  was  used  for 
all  the  derivations  presented.  Specifically,  some  of  the 
symbols  used  are: 

E  =  Electric  field  vector 
H  =  Magnetic  field  vector 
r  =  posistion  vector 
■  =  dot  product 

X  =  cross  product 
£  =  Laplacian  operator 

r  =  distance  from  z  axis 
6  =  angular  rotation  from  x  axis 

$  =  z  component  of  E  or  H  field 
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e 

o 


permittivity  of  free  space 
permeability  of  free  space 
relative  permittivity 
relative  permeability 
index  of  refraction 
conductivity 
wave  number 
radian  frequency 


Summary 

An  approximation  for  the  RCS  of  a  two  dimensional, 
perfectly  conducting  cylinder  coated  with  a  lossy  dielectric 
material  will  be  derived  using  two  equivalent  OSRC  methods. 
The  results  from  this  analysis  will  be  compared  to  the 
series  solution  and  the’  range  of  validity  of  the  OSRC 
technique  will  be  investigated.  In  addition,  the  OSRC 
solution  for  this  problem  will  be  shown  to  reduce  to  the 
OSRC  solution  for  a  perfectly  conducting  cylinder  as  the 
thickness  of  the  dielectric  layer  goes  to  zero.  The  OSRC 
technique  has  shown  excellent  results  when  applied  to 
several  two  dimensional  perfectly  conducting  bodies.  The 
OSRC  approximation  was,  in  general,  accurate  to  within  1  dB 
using  the  second  order  boundary  operator. 
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II. 


Review  of  Literature 


Introduction 

Applying  a  radiation  condition  at  infinity  has  been  used 
for  decades  to  solve  for  the  scattered  fields  from  objects 
where  the  geometry  allows  the  investigator  to  derive  an 
analytical  solution.  For  geometries  where  a  numerical 
solution  is  required  a  finite  region  must  be  artificially 
introduced  to  make  the  problem  feasible  for  computer 
application.  Over  the  years,  several  techniques  have  been 
developed  which  accomplish  this  feat  with  reasonable  degrees 
of  accuracy  and  varying  degrees  of  difficulty.  The  OSRC 
technique  consists  of  the  application  of  higher  order 
differential  operators  to  the  far  field  expression  for  the 
scattered  field  to  truncate  the  numerical  solution. 
Kriegsmann  and  others  noted  that  the  use  of  higher  order 
differential  operators  appears  to  fall  into  two 
categories.  The  first  uses  the  asymptotic  behavior  of  the 
appropriate  Greens  function  to  establish  a  series  of 
operators  that  annihilate  the  first  n  terms  of  the 
asymptotic  series  of  the  scattered  field  (Kriegsmann  and 
others,  ''987;]53).  They  describe  the  second  category  as 
follows : 

The  second  category  .  .  .  derives  an  approximate  one¬ 
way  wave  equation  in  Cartesian  coordinates  by  factoring 
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the  dispersion  relation  of  the  full  wave  equation,  and 
providing  a  rational  polynomial  interpolation  of  the 
resulting  scjuare  root  at  selected  wave  propagation 
angles.  This  results  in  reflection  free  passage  of  the 
plane  waves  propagating  at  these  angles  through  the 
lattice  truncation  plane.  The  number  of  reflection 
free  angles  and  their  values  can  be  selected  in  a 
systematic  manner  [Kriegsmann  and  others,  1987:153]. 


The  OSRC  technique  to  be  used  here  falls  into  the  first 
category.  Before  beginning  the  derivation  of  the  OSRC 
solution,  a  few  of  the  other  techniques  used  to  solve  this 
problem  will  be  reviewed. 


Summary  of  Other  Techniques 

This  problem  has  been  solved  using  a  vast  number  of 
different  techniques.  A  few  of  the  more  pertinent 
techniques  will  be  discussed.  The  first  technique  to  be 
presented  is  the  eigenfunction  expansion  solution  followed 
by  a  recently  published  geometrical  optics  solution.  Then 
several  of  the  more  recently  published  numerical 
approximations  will  be  presented.  These  techniques  include 
the  method  of  moments  (MM) ,  the  unimoment  method  and  a 
hybrid  finite  element  method. 

Eigenfunction  Expansion.  For  a  circular  cylinder  covered 


with  a  circular  dielectri-c  layer,  the  normal  boundary 
conditions  and  the  radiation  condition  at  infinity,  Tang 
used  the  eigenfunction  expansion  method  to  solve  the  system 


of  partial  differential  equations  for  the  backscattered 
field.  The  solution  consists  of  an  infinite  series  of 
Hankel  functions,  in  which  the  constants  are  ratios  of  sums 
of  Bessel  functions  and  their  derivatives,  evaluated  on  the 
conducting  and  dielectric  surfaces.  Tang  found  excellent 
agreement  between  this  solution  and  measured  values  for  the 
cases  presented  below  (Tang,  1957:628-630). 


ka  =  0.942, 

1.5 

< 

< 

10, 

e  =2.54 
r 

ka  =  2.0, 

3.188 

< 

kjb 

< 

10, 

e  =2.54 
r 

ka  =  3.5, 

5.58 

< 

kjb 

< 

10, 

e  =2.54 
r 

ka  =  0.942, 

2.308 

< 

’'2*^ 

< 

10, 

e  =6 
r 

ka  =  1.302, 

3.188 

< 

kjb 

< 

10, 

"r  =  ^ 

where, 

k  =  wave  number  in  free  space 
k2  =  wave  numb'ei'  in  the  dielectric 
a  =  radius  of  the  conductor 
b  =  radius  of  the  dielectric 

=  relative  permittivity  (Tang,  1957:628-630) 

Since  this  solution  is  a  strict  mathematical  solution  to  the 
partial  differential  equations  it  can  be  called  the  exact 
solution  for  the  circular  geometry  case. 

Geometrical  Optics.  Wang  has  used  a  Watsons 
transformation  on  the  eigenfunction  solution  to  approximate 
the  scattered  field  from  a  dielectric  coated  circular 
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cylinder  as  the  sum  of  a  geometrical  optics  or  ray  solution 
and  two  surface  waves.  Numerical  results  were  compared  to 
the  eigenfunction  solution  and  found  to  be  in  excellent 
agreement  for  ka  greater  than  2.0  and  the  thickness  of  the 
dielectric  coating  less  than  0.41  wavelengths  (Wang, 


1985:960-963) . 


Method  of  Moments.  The  method  of  moments  (MM)  has  been 
used  to  solve  the  dielectric  cylinder  problem  by  Wu  and 
Tsai.  They  formulated  the  problem  for  arbitrary  shaped 
conductor  and  dielectric  coating  as  a  system  of  coupled 
surface  integral  equations  derived  from  Maxwell's  equations, 
Green's  theorem  and  the  appropriate  boundary  conditions  for 
transverse  electric  (TE)  or  transverse  magnetic  (TM) 
polarization.  The  radiation  conditions  were  incorporated 
into  the  formulation  via  the  free  space  Green's  function. 
They  solved  the  system  of  the  integral  equations  using  flat 
pulse  expansion  functions  and  point  matching.  Several 
examples  of  dielectric  cylinders  and  dielectric  coated 
conducting  cylinders  were  given.  The  circular  geometry 
dielectric  coated  conducting  cylinder  solution  was  compared 
to  Tang's  solution.  (Wu  and  Tsai,  1977:519-522). 

For  the  lossless  dielectric  case,  Wu  and  Tsai  found 

excellent  agreement  between  the  MM  solution  and  Tang's 

eigenfunction  solution  for  the  backscattered  cross  section 

for  the  case  ka  =  0.942,  e  =2.54  and  2  <  kb  <  7 .  In 

'  r 
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addition  to  the  lossless  case,  they  computed  the  cross 
section  for  the  conductivity  (a)  of  the  dielectric  equal  to 


0.01  and  0.05  mho/meter.  It  was  concluded  that  the  effect 
of  a  lossy  dielectric  was  to  significantly  dampen  RCS 
oscillations  (Wu  and  Tsai,  1977:522). 

Another  example  that  is  relevant  here  is  the  circular 
cylinder  enclosed  in  a  dielectric  wedge  semi-circle.  Figure 
1  shows  the  geometry  and  MM  solution  for  this  problem.  The 
dotted  lines  indicate  regions  where  the  solution  was 
extrapolated  due  to  numerical  instabilities  in  the  matrix 
inversion.  Wu  and  Tsai  found  that  the  numerical 
instabilities  occurred  whenever  an  interior  resonance 
occurred  and  that  the  condition  number  was  an  excellent 
indicator  of  this  instability  (Wu  and  Tsai,  1977:520-523). 
This  seems  to  be  an  inherent  limitation  of  the  MM  technique. 

Unimoment  Method.  Chang  and  Mei  have  applied  unimoment 
methods  to  the  dielectric  cylinder  problem.  The  unimoment 
method  involves  enclosing  the  cylinder  in  an  artificial 
boundary,  usually  a  circle  for  two  dimensional  problems,  and 
enforcing  boundary  conditions  of  the  E  and  H  fields  on  this 
artificial  boundary  instead  of  the  actual  dielectric 
boundary.  The  total  field  inside  of  this  boundary  is 
expressed  as  a  finite  sum  of  unknown  constants  multiplying  a 
set  of  expansion  functions  and  the  scattered  field  outside 
of  this  boundary  is  approximated  by  a  finite  series  of 
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Figure  1.  Backscattered  RCS  of  a  Conducting  Circular 
Cylinder  Enclosed  in  a  Dielectric  Wedge-Semicircle 

(Wu  and  Tsai, 1977:522) 

Hankel  functions.  The  incident  field  is  generally  assumed 
to  be  a  plane  wave.  For  two  dimensional  problems  the 
obvious  choice  of  expansion  functions,  for  the  field  inside 
of  the  boundary,  that  satisfy  periodicity  conditions  are 
sinusoids.  Finite  element  methods  are  used  inside  of  the 
artificial  boundary  to  solve  for  the  total  field  on  the 
surface  of  the  boundary.  The  continuity  of  the  E  or  H 
fields  across  this  boundary  is  used  to  relate  the  constants 


of  the  two  series  expansions,  thus  solving  for  the  scattered 
field  or  scattering  cross  section.  Chang  and  Mei  used  this 
technique  to  solve  for  the  scattered  field  from  an  off 
center,  circular,  dielectric  cylinder  and  other  geometries 
and  found  excellent  agreement  between  their  solution  and  the 
normal  series  solution  obtained  from  separation  of  variables 
(Chang  and  Mei,  1976:35-39). 

Hybrid  Finite  Element  Methods.  As  illustrated  by  Chang 
and  Mei  above,  finite  element  methods  can  be  used  to  solve 
field  problems  in  unbounded  regions.  However,  imposing 
artificial  boundaries  is  not  efficient  when  dealing  with 
noncircular  geometries  (Jin  and  Liepa,  1988:50).  As 
implemented  by  Jin  and  Liepa,  the  hybrid  finite  element 
method  (HFEM)  encloses  tbe  scatterer  in  two  artificial 
boundaries  that  match  the  surface  of  the  scatterer.  This 
reduces  the  region  in  which  a  numerical  solution  must  be 
found.  Figure  2  shows  the  general  problem  formulation  used 
by  Jin  and  Liepa.  They  describe  the  technique  as  follov^s: 

The  basic  idea  of  the  method  is  to  draw  one  artificial 
boundary  F,  enclosing  r_  and  another  artificial 
boundary  Fg  lying  between  F^  and  F  .  To  minimize  the 
region  treated  by  FEM,  F.  is  drawn^in  such  a  way  that 
only  one  layer  of  triangular  elements  is  needed  between 
F^  and  F^.  We  will  use  FEM  in  the  region  defined  by 
F^  and  F^,  and  use  free-space  Green's  function  integral 
formulation  in  the  region  outside  F,  [Jin  and  Liepa, 
1988:51].  ^ 

Jin  and  Liepa  found  good  agreement  between  their  solution 
for  the  far  field  factor  (Bowman  and  others,  1969:6)  for  a 


Figure  2. 


Jin  and  Liepa's  HFEM  Problem  Formulation 

(Jin  and  Licpa,  1988:51) 


circular  cylinder  coated  with  a  circular  dielectric/magnetic 
material  and  the  eigenfunction  expansion  (Tang,  1957:628- 
630)  or  the  method  of  moments  with  impedance  boundary 
conditions.  In  addition,  they  compared  the  efficiency  and 
accuracy  of  the  HFEM  solution  to  a  volume  moment  method 
solution  for  the  same  geometry  but  with  a  dielectric 
material  only  (Jin  and  Liepa,  1988:52).  They  summarize  the 
results  as  follows: 

Comparison  shows  that  HFEM  has  essentially  the  same 

accuracy  as  MM  in  this  E-polarized  case.  Hov/ever, 


using  MM  in  the  case  of  H-polarization,  one  would  have 
to  compute  both  the  'fformal  and  tangential  equivalent 
electric  current  components  and  thus  have  a  larger 
resultant  matrix  .  .  .  [Jin  and  Liepa,  1988:52]. 

They  also  concluded  that  with  HFEM  the  size  of  the  matrix 

remains  the  same  if  the  coating  material  is  magnetic  as  well 

as  dielectric  while  with  MM  the  matrix  size  will  increase  by 

approximately  a  factor  of  three  (Jin  and  Liepa,  1988:52). 

Summary 

Several  different  methods  for  solving  for  the  scattering 
cross  section  of  a  dielectric  coated  conducting  cylinder 
were  reviewed.  In  general,  if  the  geometry  did  not  match  a 
set  of  orthogonal  coordinates,  a  numerical  technique  must  be 
applied  to  the  problem.  Introducing  an  artificial  boundary 
is  one  numerical  method  that  has  been  successfully  used 
(Chang  and  Mei,  1976:35-42)  and  work  has  been  done  on 
improving  the  efficiency  of  this  technique  (Jin  and  Liepa, 
1988:50-54).  However,  the  OSRC  method  takes  advantage  of 
the  natural  boundaries  of  the  system  to  derive  an 
approximation  for  the  cross  section  or  surface  currents. 

Its  advantage  is  simplicity  and  a  consistent  ease  of 
application  for  arbitrary  convex  cylinders  (Kriegsmann  and 
others,  1987:160).  Application  of  the  OSRC  method  to  a 
dielectric  coated  cylinder  is  not  evident  in  the  literature 
and  should  prove  to  be  a  natural  extension  to  the  previously 
published  work  with  this  technique. 
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III.  Problem  Formulation  and  Solution 


Introduction 

The  radar  cross  section  (RCS)  of  a  two  dimensional  object 
can  be  thought  of  as  the  'cross  section  per  unit  length  and 
has  been  called  the  scattering  width  (Knott  and  others, 
1985:157)  or  the  echo  area  (Peters,  1958:133).  For  our 
purposes,  RCS  will  mean  the  two  dimensional  RCS  and  is 
defined  as  (Bowman  and  others,  1969:7) : 


a  = 


iXi 

(  T^i  I  2 


H^l  ^ 

I  „i  I  2 


(3.1) 


where, 

r  =  radial  distance  from  the  z  axis 

s  .  . 

E  =  scattered  electric  field 

H^=  scattered  magnetic  field 
E^  =  incident  electric  field 
H  ^  =  incident  magnetic  field 


Since  the  incident  field  is  assumed  to  be  a  unit  magnitude 
plane  wave,  the  problem  of  determining  the  RCS  reduces  to 
finding  an  expression  for  the  scattered  fields. 

Starting  with  Maxwell’s  equations,  the  equations  for  the 
fields  scattered  from  a  dielectric  coated,  infinite  cylinder 
will  be  formulated  for  arbitrary  geometries  using  transverse 


magnetic  (TM)  and  transverse  electric  (TE)  polarizations. 

TM  polarization  implies  the  incident  electric  field  is 
polarized  parallel  to  the  cylinder  axis  only  while  TE 
polarization  implies  the  incident  magnetic  field  is 
polarized  parallel  to  the  cylinder  axis.  These  equations 
can  be  solved  exactly  for  a  right  circular  cylinder  coated 
with  a  circular  layer  of  dielectric  material  and  this 
solution  will  be  presented.  The  OSRC  approximation  to  the 
scattered  fields  will  be  derived  for  this  special  geometry 
and  this  approximate  solution  will  be  compared  to  the  exact 
solution.  The  cylindrical  coordinate  system  will  be  used 
throughout  this  analysis.. 

Problem  Formulation 

Figure  3  shows  the  arbitrary  geometry  problem,  a 
perfectly  conducting  cylinder  coated  by  a  layer  of 
dielectric  material  with  relative  permittivity  e^,  relative 
permeability  and  conductivity  a.  E  and  H  are  the  total 
electric  and  magnetic  fields  in  Region  I,  free  space.  e'^ 
and  are  the  electric  and  magnetic  fields  in  Region  II, 
the  dielectric  material.  The  fields  in  Region  I  consist  of 
an  incident  plane  wave  traveling  in  the  x  direction  (E^,  H^) 
and  the  fields  scattered  _from  the  dielectric  coated  cylinder 
(E  /  H  ) .  IS  the  surface  of  the  dielectric/conductor 

boundary  and  is  the  surface  of  the  dielectric/free  space 


boundary . 


# 


Assuming  an  exp(-jwt)  time  dependence,  Maxwell's 
equations  can  be  written  .as: 

curlE  =  and  curlH  =  -(a  +  jwe^e^)E 

Taking  the  curl  of  both  equations,  substituting  the  original 
equations  into  the  new  equations  and  rearranging  terms 
yields  the  homogenous  vector  wave  equation  for  the  electric 
or  magnetic  field  in  an  linear,  homogenous,  isotropic  medium 
(Harrington,  1961:37). 

curlcurlE  -  (kn^)^  E  =  0  (3.2a) 

curicurlH  -  (kn^)^  H  =  0  (3.2b) 

where, 

1/2 

k  =  wave  number  =  w(e  n  )  ' 

o  o 

.  .  1/2 
n  =  complex  index  of  refraction  =  (e  u_(  1  +  Ta/(we^e_)}  ' 
o  IT  .  ^  r  r  j  !  ^  r  o 


w  =  radian  frequency 

For  TM  polarization,  the  electric  field  is  expressed  as 
E  =  (  0,  0,  $  )  and  the  electric  field  equations  are  used  to 
solve  for  the  RCS.  For  TE  polarization  the  magnetic  field 
is  expressed  as  H  =  (  0,  0,  $  )  and  the  magnetic  field 
equations  are  used  to  solve  for  the  RCS.  Substituting  these 
field  representations  into  the  vector  wave  equations  above 
yields  the  scalar  Helmholtz  equation  for  both  polarizations. 

£$  +  (kn^)^  $  =  0  (3.3) 

where  £  is  the  Laplacian  operator  and  $  is  the  z  component 
of  the  electric  field  for  TM  polarization  or  $  is  the  z 
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component  of  the  magnetic  field  for  TE  polarization.  It  is 
implicit  in  Equation  3.3  that  divE  =  divH  =  0  (e.g.  the 
medium  has  zero  charge  density) . 

In  Region  II  $  =  where  is  the  total  field  in  the 
dielectric  material  and  satisfies  Equation  3.3.  In 
Region  I,  $  =  and  n^  =  1,  where  and  are  the 

scalar  components  of  the  incident  and  scattered  fields 


respectively.  Substituting  into  Equation  3.3  yields: 


=  0 


(3.4) 


However,  since  exists  whether  or  not  the  scattering 
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object  is  present,  it  must  satisfy  Equation  3.3  by  itself. 
Therefore,  Equation  3.4  reduces  to: 

=  0  (3.5) 

In  addition  to  Equations  3.3  and  3.5,  the  boundary 
conditions  at  and  are  needed  to  find  a  unique  solution 


for  the  scattered  fields. 

The  boundary  conditions  at  are: 

n  X  E  =  0  (3.6a) 

n  X  H  =  J®  (3.6b) 

n  •  E  =  Pg/e^  (3.6c) 

n  •  H  =  0  (3.6d) 

The  boundary  conditions  at  are: 

n  X  E  =  n  X  (E^  +  E^)  =  n  X  E^  (3.7a) 

n  X  H  =  n  X  (H^  +  H^)  =  n  X  (3.7b) 

n  ■  E  =  n  ■  (E^  +  E®)  =  Cj.  (n  •  E^)  (3.7c) 

n  ■  H  =  n  ■  (H^  +  H®)  =  |ij.(n  •  H^)  (3.7d) 

where. 


n  =  unit  vector  normal  to  the  respective  surface 

g 

J  =  surface  current  vector 
Pg  =  surface  electric  charge  density 
For  incident  TM  polarization.  Equations  3.6a,  3.7a  and  3.7b 
reduce  to: 

=  0  on  (3.8a) 

on  (3.8b) 
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(3.8c) 


on  r. 


For  incident  TE  polarization,  equations  3.6b,  3.7a  and  3.7b 
reduce  to: 


on  r. 


(3.9a) 


on  r. 


(3.9b) 


on  r. 


(3.9c) 


Equations  3.3,  with  i  3.5  and  3.8  form  the  complete 

problem  for  incident  TM  polarization,  while  Equations  3.3, 
3.5  and  3.9  form  the  complete  problem  for  incident  TE 
polarization. 

In  order  to  implement  the  OSRC  technique  it  is  necessary 
to  express  Equation  3.5  in  its  integral  equation  form  using 
Green's  Theorem  and  the  two  dimensional  free  space  Green's 
function  for  the  scalar  Helmholtz  equation.  Using  these 
techniques,  Equation  3.5  can  be  recast  as: 


[  G(r,r') 


_  $  s  SG{r,r'_)  ^  p  ^g) 


where , 


G(r,r')  =  Green's  function  = 


(k|r  -  r'  1) 


and  H'  MX)  is  a  Bessel  function  of  the  third  kind  or 
n  '  ' 

Hankel  function.  The  boundary  conditions  (Equations  3.8b 
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and  3.8c  for  incident  TM  polarization  or  Equations  3.9b  and 
3.9c  for  incident  TE  polarization)  can  be  used  to  substitute 
for  the  scattered  field  and  its  normal  derivative  in  the 
integrand  of  Equation  3.10.  Thus,  the  expression  for  the 
scattered  field  becomes: 

.D 


[  G(r,r')  I  ^ 
^r  5n 


jD  iG [rj rl).  j 

5n 


[  G(r,r') 


Sn 


i  6Gfr,r' ) 
(Sn 


]  6s'  (3.11) 


where  p  is  u  for  TM  polarization  or  e  for  TE 
polarization. 


Exact  Solution  for  Circular  Geometries 

For  a  circle  of  radius  b  and  a  circle  of  radius  a, 

where  b  >  a,  the  normal  derivatives  in  all  of  the  above 
equations  become  derivatives  with  respect  to  r' ,  the  radial 

coordinate.  Thus,  Equation  3.11  becomes: 

,27r 


=  b 


[  G(r,r') 


6$ 


D 


r  6r' 


6Gir,r'.l  ^ 
6r' 


n27r 


-  b 


[  G(r,r') 


6±Z 

Sr’ 


$ 


i  SG(r,r' ) 
Sr' 


]  60'  (3.12) 


Since  we  are  interested  in  the  behavior  of  the  scattered 
field  as  r  goes  to  infinity,  the  first  term  of  the 
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asymptotic  expansion  of  the  Green's  function  and  its 
derivative  can  be  used  in  Equation  3.12.  These  expansions 
are: 


G(r,r') 


-jkr'cos(e  -  0') 


^G(r,rM 

Sr 


-1/-  -jkr'cos(e  -  e') 

-3kGQCOs(e  -  O')  e  ' 


where , 


j  (kr  +  n/4,) 

172 


(3.12a) 

(3.12b) 

(3 . 12c) 


(STrkr) 

Substituting  Equations  3.12a  and  3.12b  into  Equation  3.11 
yields: 


27r 


D 


r  1  .D,  ~jkbcos(e  -  ©') 

[  —  —  +  jkcos(0  -  0')  $  ]  e  ^  ' 

Pr  Sr' 


de‘ 


—  +  jkcos(0 

Sr' 


0') 


g-*jkbcos(0  -  0')  ^0, 
(3.13) 


s  s 

where  is  the  far  field  approximation  for  $  . 

The  incident  field  is  assumed  to  be  a  unit  magnitude 
plane  wave  traveling  in  the  positive  x  direction.  This  can 
be  expressed  as  the  normal  complex  exponential  function  or 
as  a  series  of  cylindrical  wave  functions  (Harrington, 


1961:231) ; 


.  .  . 

.1  jkx  lkrcos(0)  „  .n  ^  ,  in0 

§=e-'=e'’  (kr)  e-* 

n=-oo 

Now  it  is  necessary  to  solve  Equation  3.3  to  find  an 


(3.14) 
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expression  for  to  use  in  the  integrand  of  Equation  3.13. 

The  method  of  eigenfunction  expansion  will  be  used  to 
D  s 

solve  for  $  and  $  for  either  polarization.  Any  linear 
combination  of  Bessel  functions  will  suffice  as  an 
eigenfunction  for  Equation  3.3  (Baker,  1988) .  However, 
Hankel  functions  were  chosen  because  they  exhibit  the 
correct  behavior  for  traveling  wave  solutions.  Thus,  the 
eigenfunctions  are: 

(x)  =  outward  going  wave  solution  =  +  j 

(2) 

^  (x)  =  inward  going  wave  solution  =  **  ^ 

where  J  (x)  are  Bessel  functions  of  the  first  kind  and  Y  (x) 
n  n 

are  Bessel  functions  of  the  second  kind. 

In  additional  to  the  boundary  conditions,  all  the  field 
solutions  must  satisfy  periodicity  conditions  with  respect 
to  and  the  scattered  field  must  satisfy  the  Sommerfeld 
radiation  conditions. 

Thus  the  solution  for  may  be  written  as: 

n=oo 

=  S  j”  (  A  H^^^(kn  r)  +  B  (kn  r)  )  e^^®  (3.15) 

„','nn'o'  nn'o'^  '  * 

n=-oo 

and  the  solution  for  may  be  written  as: 

\  j"  C  Hj(^)(kr)  e^""®  (3.16) 

n=-oo 

where  A  ,  B  ,  and  C  are  the  unknown  coefficients, 
n  n  n 
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Enforcing  the  boundary  conditions  at  and  yields: 


B  =  -  A 

n  n 


B  =  -  A 

n  n 


H<2> (kn^a) 


(kn^^a) 


for  TM  polarization 


for  TE  polarization 


(3 . 17a) 


(3.17b) 


(1) 


,(1) 


H^!;"^(kb)  J^(kb)  -  H;^"^(kb)  J^(kb) 

"  n  (kb)  -  p  (kb)  ^ 


r  n  n 


(3 . 18) 


C  = 
n 


Pr°n  '  "o 

"c  “n  "a"’  -  pA  <>*> 


(3.19) 


where,  F(z)  implies  the  derivative  of  F(z)  with  respect  to 
the  argument.  For  TM  polarization: 

=  Hj^^^knQb)  H^2^(knQa)  -  (kn^b)  H^^^^kn^a)  (3.20a) 

=  H^^^kn^b)  H^^^kn^a)  -  (kn^b)  (kn^a)  (3.20b) 


n  n 


n 


n 


n 


H<2> (kn^a) 


While,  for  TE  polarization: 


n 


H<=> (kn„a) 


(3.20c) 


(kn  b)  (kn  a)  -  (kn  b)  (kn^a)  (3.21a) 

n  n'o'n'o'  n'o'n'o'  ' 

(kn  b)  (kn  a)  -  (kn  b)  (kn  a)  (3.21b) 

n  n'o'n'o'  n'o'n'o''  ' 


(3.21c) 
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It  is  possible  to  generate  at  least  two  equivalent 
expressions  for  the  RCS  based  on  the  equations  presented 
above.  Having  solved  for  the  coefficients  in  the  series 
expansions,  Equation  3.16  can  be  substituted  into  Equation 
3.1  directly  or  Equation  3.15  can  be  substituted  into 
Equation  3.13  and  the  resulting  integral  expressions 
evaluated  before  substituting  into  Equation  3.1.  The  latter 
method  will  be  used,  even  though  it  requires  more  work, 
because  it  generates  a  series  expression  for  the  RCS  which 
is  easily  compared  to  an  OSRC  approximation. 

Substituting  the  expressions  for  and  into  Equation 
3.13  and  evaluating  the  resulting  integrals  yields  the 
following  expression  for  the  RCS. 


(irkh) 

^ 


2  1  . 

^  “n  ^n^  ^o  K  ‘  ^r  ^n  )cos(n0) 

n=0 


(3.22a) 


where,  is  1  for  n  =  0  and  2  otherwise  and. 


and  p^. 


T?  — 

Jj,(kb)  - 

j„(kb) 

(3 . 22b) 

£i  — 

n 

p  n 

^r  o  n 

-  d" 

'^n  n 

II 

W  P 

H^^^(kb)  / 

(kb) 

(3.22c) 

1  2 

D^,  and  are  polarization  dependent  and  are 


defined  above. 


OSRC  Approximation  for  Circular  Geometries 

Two  equivalent  OSRC  solutions  will  be  derived.  The  first 
method  produces  a  series  expression  for  the  RCS  that  is 
easily  compared  to  the  exact  solution  presented  above.  The 
other  method  shows  how  the  OSRC  solution  for  the  dielectric 
coated  cylinder  is  an  extension  of  the  OSRC  solution  for  the 
perfectly  conducting  cylinder  and  is  more  amenable  to 
numerical  solutions  for  arbitrary  geometries.  All  of  the 
work  presented  here  uses  the  second  order  boundary  value 
operator  used  by  Kriegsmann  and  others  in  their  development 
of  the  OSRC  approximation  for  the  RCS  of  a  perfectly 
conducting  cylinder  (Kriegsmann  and  others,  1987:155). 

For  a  circle  of  radius  a,  a  circle  of  radius  b, 
where  b  >  a,  and  using  the  far  field  approximations  for  the 
Green's  function  and  its  normal  derivative  presented  above, 
Equation  3.10  can  be  reduced  to: 

27r 


*£  =  t-G„ 


[  ^  +  jkcos(e  -  S')  -  ®')d9' 

Sr’ 


(3.23) 


s  • 

where  is  the  far  field  approximation  for  the  scattered 
field  and  is  defined  in  Equation  3.12c. 

The  first  OSRC  method  uses  the  boundary  conditions  to 
substitute  for  the  scattered  field  and  its  derivative  in  the 
integrand  of  Equation  3.23  in  terms  of  the  incident  field 
and  the  field  in  the  dielectric  material.  Then  the  OSRC 
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♦ 


operator  is  used,  instead  of  the  Sommerfeld  radiation 
conditions,  to  find  an  approximation  for  the  field  in  the 
dielectric  material.  The  resulting  integrals  are  evaluated 
and  the  final  expression  substituted  into  Equation  3.1. 

This  method  is  very  similar  to  the  procedure  presented  above 
for  the  exact  solution  and  produces  a  series  expression  for 
the  RCS  which  is  almost  identical  to  Equation  3.22. 

The  second  method  uses  the  OSRC  operator  to  find  an 
approximate  expression  for  the  normal  derivative  of  the 
scattered  field  in  the  integrand  of  Equation  3.23  and  for 
the  field  in  the  dielectric  material.  The  boundary 
conditions  are  used  to  substitute  for  the  scattered  field  on 
the  surface  of  the  dielectric  material  and  the  resulting 
integrals  are  evaluated'  to  find  the  far  field  approximation 
for  the  scattered  field.  The  expression  for  the  RCS  is 
found  by  substituting  into  Equation  3.1.  This  method 
produces  a  solution  which,  for  the  special  case  of  a  =  b, 
can  be  easily  shown  to  reduce  to  the  OSRC  approximatio'.  for 
the  RCS  of  a  perfectly  conducting  cylinder  published  by 
Kriegsmann  and  others  (Kriegsraann  and  others,  1987:156). 
Since  both  methods  use  the  OSRC  approximation  for  the  field 
in  the  dielectric  material,  the  details  of  this  derivation 
will  be  presented  first. 

The  assumed  form  for  and  the  boundary  condition  on 
remain  the  same  as  for  the  exact  solution.  Therefore, 
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can  be  expressed  as  a  summation  of  eigenfunctions  with  one 
unknown  coefficient. 
n=oo 

=  L  A  [  (kn  r)  -  B  (kn  r)  ]  e^^®  (3.24a) 

in'-n'o'  nn'o'-*  '  ' 

n=— 00 

where,  •  - 

(kn  a) 

B  =  - ryr -  for  TM  polarization  (3.24b) 

Hf>(kn„a) 

or, 

H  ^  (kn  a) 

B  =  — - ^  for  TE  polarization  (3.24c) 

"  ^ (kn  a) 

n  '  o  ' 

and  A  is  the  unknown  coefficient.  A  is  found  by  using  an 
n  n 

n^^  order  OSRC  boundary  operator  and  the  boundary  conditions 
on  T^. 

The  n^^  order  OSRC  boundary  operator,  is  used  in  lieu 

of  the  Sommerfeld  radiatron  condition  by  imposing  the 
following  condition  on  the  solution: 

B^{  =  0  on  (3.25) 

Equation  3.25  is  an  approximation  with  an  error  term  that  is 

proportional  to  r^  i/2)^  where  r  is  the  distance  from  the 
origin  to  the  surface  of  the  dielectric  material  and  n  is 
the  order  of  the  operator.  Appendix  A  explains  the 
motivation  for  using  this  approximation  and  presents  a 
derivation  of  the  error  term  and  the  second  order  operator. 
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The  second  order  OSRC  op&rator  is  (Kriegsmann,  1987:155); 


B, 


-  jk  - 


60' 


'2  6  r  '  2r  2r  (1  -  jkr) 

Substituting  Equation  3.26  into  Equation  3.25  and 
rearranging  terms  yields: 

,s 


6$: 

Sr 


k  B  {  $  }  on  r. 


where. 


B 


3  - 


60' 


_1 

4 


2kr 


2kr(l  -  jkr) 


(3.26) 


(3.27a) 


(3.27b) 


Rearranging  the  boundary  conditions  by  solving  Equations 
3.8b  and  3.8c  or  Equations  3,9b  and  3.9c  for  the  scattered 
field  and  its  derivative  and  substituting  the  resulting 
expressions  into  Equation  3.27a  yields: 


1_ 

P, 


6$ 


D 


5r 


-  kB  {  $^  )  = 


6§' 


6r 


-  kB  {  4>  } 


(3.28) 


where  p^  is  for  incident  TM  polarization  or  for 


D 


TE  polarization.  Substituting  the  assumed  form  for  $  and 
multiplying  both  sides  of  Equation  3.28  by  exp(-jm0)  and 
integrating  from  0  to  27r  yields  the  unknown  coefficient 


A  = 
n 


J^(kb)  -  J(kb) 

n  '  p 

on  ^r  '^n  n 


P  P 
'^r  n 


(3.29a) 


1  2 

where  D^,  p^  and  are  polarization  dependent  and 
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defined  above  in  Equations  3.20  and  3.21  and: 


2 

1  n  -  1/4  ooKN 

^n  ^  “  2kr  "  2kr(l  -  jkr)  (3.29b) 

Now  that  the  OSRC  approximation  for  the  field  in  the 

dielectric  layer  is  known,  the  approximations  for  the 

scattered  field  can  be  found  using  the  two  methods  described 

above. 

Using  the  boundary  conditions  to  substitute  for  the 
scattered  field  and  its  derivative  in  the  integrand  of 
Equation  3.23  yields: 


"?=  bGo 


,27r 


[1  ^  jKoosO  -  6')  #°  J  -  ®'>d9' 


r  ir' 


p2ir 


-  bG, 


t  ili  +  jkcos(e  -  9') 

5r' 

0  (3.30) 


which  is  identical  to  Equation  3.13  above.  The  two  terms  in 
the  second  integral  of  Equation  3.30  cancel;  therefore,  its 
contribution  is  zero.  Substituting  the  expression  for  in 
the  first  integral  yields: 

n=<» 


=  kbGo  S  j"  I 
n=-<» 


(3.31a) 


where, 

.27r 


I  = 


r  r>2  ,  .  r^l  rvfVT  3  n0  '  -  3kbCOS(0  -  0  ).o, 

[n  D  +  j  p  D  cos(0  -  0  )1  e-*  '  'd0' 

I  (3.31b) 
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The  integral  in  equation  *3.31  is  evaluated  using  the 
following  Bessel  function  identities. 

r2n 


J^(z)  - 


and 


j  2.j-"  e^^^®  j^(z)  = 


gjne'->  jzcos(0  -  0') 


(3.32) 


-27r 


cos(0  -  0')  e 


jn0'-  jzcos(0  -  0') 


d0' 


0  (3.33) 

Substituting  Equations  3.32  and  3.33  into  Equation  3.31 
yields  the  OSRC  approximation  for  the  scattered  field. 
n=oo 

=  27rkbG  S  A  [  n  J  (kb)  -  p  dJ  J  (kb)  ]  e 

n=-oo  .  ”  r  n  n  (3.34) 

Substituting  into  Equation  3.1  yields  the  OSRC  approximation 
for  the  RCS  using  the  first  method. 

2 


_  (TTkb) 


n=oo 


Z  a  A„(  n  J  (kb)  -  p  J  (kb)  )cos(n0) 
^^nn'  onn'  '  ^rnn'  '  '  '  ' 

n=0 


(3.35) 


where  is  1  for  n  =  0  and  2  otherwise.  A  thorough 


comparison  of  this  approximation  and  the  exact  solution  is 
presented  below. 

The  second  method  uses  the  OSRC  boundary  operator  and  the 
t  •  s  D  i 

boundary  condition  $  =  i  -  $  on  the  surface  of  the 

dielectric  layer  regardless  of  polarization  to  find  a 

% 

substitution  for  the  derivative  of  the  scattered  field  in 
the  integrand  of  Equation  3.23. 
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(3.36) 


Si  *  s  "  D  i 

=  k  B  (  $  }  =  k  B  {  $  -  #  }  on  r2 

Substituting  Equation  3.36  and  the  boundary  condition 
mentioned  above  into  Equation  3.23  yields: 


$f  =  kbG^  (  -  I2  ) 


(3.37a) 


where, 


p27r 

=  [  B{  +  jcos(0  -  e')  “  ®'^d0' 

J  0  (3 . 


(3 . 37b) 


I2  =  [  B(  +  jcos(0  -  0')  ]  e 


1  ^-jkbcos(0  -  e') 


(3.37c) 


Using  the  series  expression  for  and  the  exponential 
expression  for  in  the  argument  of  the  B  operator  yields: 


B(  =  2  j”  A  e^^®  on  F, 

^  ^  ^  ^  ‘ 
n=-oo 


(3.38a) 


^  ^  2kb  8kb(l^-  jkb) 


jkbcosG  +  (kbsinO)  .  jkbcos0 
2kb(l  -  jkb)  > 


(3 . 38b) 


Using  Equation  3.38a, 


1=2  j^A  D^(  [/3®  +  jcos(0  -  0')]  jkbcos(0  -  0')^q, 

^  Jq  (3.39) 

Which  is  easily  evaluated  using  Equations  3.32  and  3.32  to 


give: 
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n=oo 

=  27r  S  A^D^(/3°J^(kb)  -  J^(3cb))e^^® 


n=-oo 


n=oo  .  _ 

•2-k  2  p  E 
n=-co  ^  ^ 

(3.40a) 

(3.40b) 


where, 


[  Jn(kb)  - 


n 


n  -  p 

on  n 


Evaluating  required  a  little  more  work.  Substituting 
Equation  3.38b  into  3.37c  yields: 

.27r 


X  r  •  1  1 

I-=  [:  -  ^TT-T-  + 


2kb  8kb(l-jkb) 


] 


jkb(cose'  -  cos(©  -  ©')l!,c,, 


jkb 


2kb(l  -  jkb) 


27r 


COS0'  e 


»  ^jkb(cos0'  -  cos(0  -  e'))^Q, 


iiskl 


4kb(l  -  jkb) 


27r 


(1  -  cos(20'))  -  cos(G  -  0')) 


d0' 


-27r 


+  j 


cos(0  -  0')  "  ®'^^d0'  (3.41) 


.._2. 


Where  the  identity  sin  0  =  (  1  -  cos (20)  )/2  has  been  used. 
Using  the  results  of  a  good  integral  table  (Gradshetyn  and 
Ryzhik,  1980:482)  and  combining  terms,  Equation  3.41  reduces 
to: 

I2  =  27r  [  Kq  Jq(z)  +  sin(0/2)  J^(z)  +  K2  cos0 

(3.42a) 

where  z  =  2kbsin(0/2)  and 
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K  =  T  -  -1-  _  kb 

0  2kb  4(1  -  jkb) 


8kb(l-jkb) 


(3.42b) 


IT  —  1  _ ± _ 

1  2(1  -  jkb) 


kb 


^2  4(1  -  jkb) 


(3.42c) 


(3 .42d) 


Using  Equations  3.40  and  3.42  the  final  result  for  the  OSRC 
approximation  for  the  scattered  field  using  the  second 
method  is: 

=  -27rkbGQ  [  Kq  Jq(z)  +  sin(e/2)  J^(z)  +  cos0  J2(z) 


+ 


n=«)  ,  ^ 

D  S  E  1 

^r  n  n  -  ■■ 

n=-oo 


(3.43) 


Substituting  equation  3.43  into  Equation  3.1  yields  the 
final  expression  for  the  RCS: 


:  -  (jTkb) 

^2  k 


K 


Q  Jq(z)  +  sin(e/2)  Jj^(z)  +  COS0  J2(2) 


n=oo 

+  p^  S  E„  O':  cos(n0) 

""  n=0  ^ 


(3.44) 


In  order  to  show  that  Equation  3.44  reduces  to  the  OSRC 
approximation  for  the  RCS  of  a  conducting  cylinder  for  the 
special  case  of  a  =  b  and  incident  TM  polarization,  it  is 
necessary  to  show  that  the  summation  goes  to  zero  and  to 
make  a  high  frequency  approximation  for  the  terms.  For 
a  =  b  and  incident  TM  polarization,  =  0  for  all  n  and 
every  term  in  the  summation  in  Equation  3.44  is  identically 
zero.  The  high  frequency  approximation  for  the  terms  is: 
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(1-jkb) 


kb 


and  ignoring  terms  that  are  proportional  to  (Jcb) 
expressions  for  the  terms  reduce  to: 


-2 


^0  = 


3  - 


2kb 


_  1 
4 


The 


(3.45a) 


= 


^2  = 


1 

4 


2kb 


Using  the  Bessel  function  identity: 


J2(Z) 


(2/z)J^(z)  -  Jq(z) 


(3.45b) 


(3.45c) 


and  the  trigonometry  identity  1  +  cos0  =  2  cos  0,  the 
expression  for  the  RCS  reduces  to: 


*  _  (^b), 

^3  k 


[  j  - 


_  1 


2kb  ■  2  c°s‘(0)  ]  Jq(z)  +  G  J^(z) 

(3.46a) 

G  =  [1  +  j/(2kb)]  sin(0/2)  +  cos(0)/2z  (3.46b) 

Equation  3.46  is  identical  to  the  OSRC  approximation  for 
the  RCS  of  a  conducting  cylinder  published  by  Kriegsmann  and 
others  (Kriegsmann  and  others,  1987:156). 


Comparison  of  the  OSRC  and  Exact  Solutions 

Two  expressions  comparing  the  OSRC  solution  using  the 
first  method  and  the  exact  solution  were  derived.  The  first 
expression  shows  how  the  OSRC  solution  converges  to  the 
exact  solution  as  kb  gets  large.  The  second  expression  shows 
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how  the  properties  of  the  dielectric  layer  play  a  large  role 
in  determining  the  overall  accuracy  of  the  OSRC  solution. 

Comparing  Equation  3.22  to  Equation  3.35,  the  only 
difference  occurs  for  the  /3  terms  in  the  coefficients  E 


n 

and  A^.  Repeating  the  definition  of  these  terms: 


n 


3  - 


n  2kr 

for  the  OSRC  solution  and, 

oE 


n  -  1/4 
2kr(l  -  jkr) 


n 


(kb)  /  (kb) 


(3.29b) 


(3.23b) 


for  the  exact  solution.  Using  the  asymptotic  expansions  for 
the  Hankel  function  and  its  derivative  (Abromowitz  and 
Stegun,  1972:364-365)  the  first  four  terms  in  the  ratio  of 
Equation  3.23b  are: 


(kb) 

(kb) 


3  - 


2kb 


+  3 


,  1  -  4n 


8(kb) 


1  -  4n 
8(kb)' 


+  0(kb) 


-4 


(3.47) 


Rewriting  the  last  term  in  Equation  3.23b  as  a  converging 
geometric  series: 

2 


oO  •  1  .  1  -  4n  ,  .  .  1 

8(kb>2  ^ 


kb 


(kb)^  (kb)^ 


+  .  .  .  ] 


(3.48) 

Comparison  of  Equations  3.47  and  3.48  shows  that  the  OSRC 
approximation  should  converge  to  the  exact  solution  as  kb 
gets  large.  This  convergence  should  occur  fairly  quickly 
since  the  two  series  agree  to  at  least  order  three. 
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However,  the  overall  accuracy  of  the  OSRC  technique  is 
determined  by  the  difference  between  each  coefficient  in  the 


two  series  solutions.  Letting  the  exact  solution  be 
represented  as: 

2 


_  (7rkb) 


n=oo 

S  a  F  cos(n0) 
n=0 


(3.49) 


and  the  OSRC  solution  as: 

2 


:  ^  (TTkb) 

k 


n=«>  Q 

2  a  F  cos(n0) 
n=0 


(3.50) 


then  the  accuracy  of  the  OSRC  technique  is  determined  by  how 
0  E 

closely  F^  approximates  F^  for  all  n  in  which  the  magnitude 
E 

of  Fj^  is  significant.  Jn.  order  to  derive  an  expression 

O  E 

which  relates  F  to  F  ,  let: 

n  n 


=  e  and 
'^n  '^n  n 


Then  F  can  be  written  as: 
n 


F^  = 

n  n  n 


[n^JnCkb)  -  PrDX(to)]  "  2  "  "e  1  .  ^  ' 


n  D  -  p  J0  D  +  p  e  D 
o  n  *^r'^n  n  n  n 


(3.51) 


Rearranging  terms,  Equation  3.51  becomes: 


F^  = 
n 


n 


[n  D^J  (kb)  -  p  D^J  (kb)  ]  e  J  (kb) 
••onn'  '  ^rnn'  nn'  ' 


(3.52) 


+ 


p  e  D 
^r  n  n 
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After  some  manipulation,  equation  3.52  can  be  rearranged  as: 


F‘^  =  F^{l  +  e(l  +  Y)  (1-X  +  X^- 

n  n  '  n'  n^  '  n  n 


.  .)}  (3.53a) 


where. 


Y  = 
n 


[  n  D^J  (kb)  -  p  D^j  (kb)  ]  p  J  (kb)  -  p  D^F® 


(3.53b) 


X 


e  p  D 
n  ^r  n 


n 


n  -  p  p^ 
o  n  '^r  n  n 


(3 . 53c) 


Due  to  the  complexity  of  the  Y^^  and  X^  terms  in  Equation 
3.53,  it  is  extremely  difficult  to  predict  how  closely  F*^ 

E 

approximates  F  .  However,  if  both  Y  and  X  are  much  less 
n  n  n 

than  one,  then  the  OSRC  technique  should  produce  a  good 
approximation  for  the  RCS  for  even  small  values  of  kb.  The 
dielectric  material  parameters  e^,  u^,  o  and  the  thickness 
of  the  dielectric  layer  should  have  a  significant  effect  on 
the  accuracy  of  the  OSRC  solution. 


Summary 

The  RCS  of  a  circular  cylinder  coated  with  a  lossy 
dielectric  and/or  magnetic  material  was  derived  using  the 
method  of  eigenfunction  expansion  and  the  OSRC  approximation 
technique  for  both  incident  polarizations.  Two  equivalent 
OSRC  solutions  were  derived.  One  was  shown  to  reduce  to  the 


OSRC  approximation  for  the  RCS  of  a  conducting  cylinder 
alone  as  the  thickness  of  the  dielectric  material  went  to 
zero.  The  other  OSRC  solution  was  compared  to  the 
eigenfunction  solution.  The  OSRC  solution  and  the 
eigenfunction  solution  were  nearly  identical;  however,  the 
overall  accuracy  of  the  OSRC  technique  was  shown  to  be 
goverened  by  a  complicated  expression  which  was  dependent 
the  dielectric  material  parameters  as  well  as  the  kb 
product. 


IV,  Results 


Introduction 

Nine  computer  programs  were  written  to  compare  the  OSRC 
solutions  to  the  eigenfunction  solution  for  three  different 
circumstances.  One  set  of  programs  computes  the  monostatic 
RCS  for  fixed  frequency  and  conductor  radius  and  increasing 
dielectric  coating  radius.  Another  set  of  programs  computes 
the  monostatic  RCS  for  fixed  geometry  and  variable 
frequency.  The  last  set  of  programs  computes  the  bistatic 
RCS  for  fixed  geometry  and  frequency  for  0  <  0  <  tt  radians. 
Each  set  of  programs  consists  of  one  program  to  compute  the 
eigenfunction  solution  and  two  programs  to  compute  the  two 
OSRC  solutions  presented  in  Chapter  3 .  Both  OSRC  programs 
produced  identical  results,  within  the  limits  of  the  graphs, 
for  all  of  the  data  presented  below.  Therefore,  only  one 
OSRC  plot  appears  on  each  of  the  figures.  A  listing  of  the 
programs  is  provided  in  Appendix  C.  All  of  the  programs  use 
a  common  set  of  subroutines  to  evaluate  the  necessary  Bessel 
functions.  A  complete  description  of  the  algorithms  and  a 
listing  of  the  subroutines  is  provided  in  Appendix  B.  For 
all  of  the  programs  the  summation  loops  were  stopped 
whenever  the  magnitude  of  the  coefficient  just  computed  was 
five  orders  of  magnitude  less  than  the  magnitude  of  the  sum. 
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The  programs  accept  complex  valued  and  /Xj.  and  a  non  zero 
conductivity  a. 

In  order  to  limit  the  amount  of  data  presented,  a 
few  cases  were  chosen  that  were  representative  of  the 
general  trends  observed.  For  the  fixed  frequency  and 
variable  thickness  programs,  a  conductor  radius  of  0.25 
wavelengths  and  1  wavelength  for  k  =  200  (ka  =  1.57  or  ka  = 
6.28)  was  chosen  as  the  baseline  for  comparison  for  both 
polarizations.  The  values  chosen  for  the  real  part  of  the 
dielectric  material  parameters,  and  /x^,  remained  the 
same  at  2.54  and  1.0,  respectively,  for  all  the  data.  The 
values  chosen  for  the  imaginary  part  of  the  relative 
permittivity  were  0.0,  0.25,  1.0  and  in  some  cases  1.5.  The 
imaginary  part  of  the  relative  permeability  was  left  at  0.0. 
For  the  fixed  geometry  and  variable  frequency  programs,  the 
kb  product  was  varied  from  0.06  to  7.1  for  b  =  3a  and  from 
0.17  to  20.5  for  b  =  8.5a  for  both  polarizations.  The 
dielectric  material  parameters  were  varied  as  in  the 
variable  thickness  case  described  above.  In  addition  to  the 
variable  thickness  and  variable  frequency  cases,  several 
representative  sets  of  input  parameters  were  chosen  for  the 
bistatic  RCS  computations.  The  conductivity  was  set  to  zero 
for  all  the  data. 

A  non-zero  imaginary  part  of  either  the  relative 
permittivity  or  permeability  implies  the  material  absorbs 


power  from  the  field.  For  this  reason,  the  case  where  the 
imaginary  part  of  was  zero  was  called  the  lossless  case; 
0.25  was  called  the  slightly  lossy  case;  1.0  was  called  the 
lossy  case  and  1.5  was  called  the  very  lossy  case. 

The  above  parameters  were  chosen  because  they  illustrate 
how  the  OSRC  solution  converges  to  the  eigenfunction 
solution  as  the  kb  product  or  the  material  parameters  vary. 
The  goal  here  is  to  determine  the  ranges  of  validity  of  the 
OSRC  technique,  not  to  model  a  specific  dielectric  coating. 
The  OSRC  approximation  was  said  to  agree  well  with  the  exact 
solution  whenever  the  difference  between  the  two  solutions 
was  less  than  2  db. 

First  the  variable  thickness  data  will  be  presented, 
followed  by  the  variable  frequency  data  and  then  the 
bistatic  data.  A  brief  summary  of  the  pertinent  features  of 
each  data  set  is  included  with  the  data.  A  complete 
analysis  of  all  the  data  is  included  following  the  bistatic 
data. 

Variable  Thickness  Data 

Figures  4  through  6  are  plots  of  the  RCS  versus  kb  for 
the  lossless,  slightly  lossy  and  lossy  cases  for  ka  =  1.57 
and  incident  TM  polarization.  Figures  7  through  9  are  the 
same  data  for  incident  TE  polarization  and  Figure  10  is  the 
very  lossy  data.  For  all  of  these  cases,  the  OSRC  solution 


Figure  4 


RCS  for  ka  =  1.57,  Lossless  Layer 


TM  Polarization,  Variable  Thickness 


Figure  5.  RCS  for  ka  =  1.57,  Slightly  Lossy 

Layer 
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Figure  10.  RCS  for  ka  =  1.57  Very  Lossy  Layer 


agreed  well  with  the  eigenfunction  (exact)  solution  for  kb 
less  than  approximately  3.0.  This  is  equivalent  to  the 
thickness  of  the  dielectric  layer  less  than  one-half 
wavelength.  When  the  thickness  exceeded  one-half 
wavelength,  the  OSRC  solution  agreed  with  the  exact  solution 
for  certain  combinations  of  loss  and  kb  only.  For 
incidentTM  polarization'  alid  the  imaginary  part  of  the 
relative  permittivity  equal  to  0.25,  Figure  5,  the  kb 
product  had  to  be  greater  than  7.2  for  the  two  solutions  to 
agree.  For  the  lossy  case.  Figure  6,  the  two  solutions 
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agreed  for  all  values  of  kb.  For  incident  TE  polarization 
and  the  imaginary  part  of  the  relative  permittivity  equal  to 
0.25,  Figure  8,  kb  had  to  be  greater  than  8  for  the  two 
solutions  to  agree.  For  the  lossy  case.  Figure  9,  the  kb 
had  to  be  greater  than  7.0  for  the  two  solutions  to  agree. 
For  the  very  lossy  case.  Figure  10,  the  two  solutions  agreed 
for  all  values  of  kb. 

Figures  11  and  12  are  plots  of  the  RCS  versus  kb  for  the 
lossless  and  slightly  lossy  cases  and  ka  =  6.28  for  incident 
TM  polarization.  For  the  lossless  case.  Figure  11,  the  OSRC 
solution  agreed  well  with  the  exact  solution  for  kb  less 
than  approximately  7.5.  This  is  equivalent  to  the  thickness 
of  the  dielectric  layer  less  than  0.2  wavelengths.  For  the 
kb  product  greater  than  14,  the  two  solutions  agreed 
everywhere  except  at  nulls.  For  the  slightly  lossy  case. 
Figure  12,  the  OSRC  solution  was  within  2  db  of  the  exact 
solution  for  all  values  of  kb. 

Figures  13  through  15  are  the  data  for  incident  TE 
polarization.  For  the  lossless  dielectric  case.  Figure  13, 
there  was  virtually  no  region  where  the  OSRC  solution  agreed 
well  with  the  eigenfunction  solution.  For  the  slightly 
lossy  case.  Figure  14,  The  two  solutions  agreed  whenever  kb 
was  greater  than  8  and  for  the  ossy  case,  Figure  15,  the 
two  solutions  agreed  for  all  values  of  kb. 
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Figure  11.  RCS  for  ka  =  6.28,  Lossless  Layer 


TM  Polarization,  Variable  Thickness 


Figure  12.  RCS  for  ka  =  6.28  Slightly  Lossy 

Layer 
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TE  Polarization.  Variable  Thickness 


a  Exact  ■¥  OSRC 


Figure  15.  RCS  for  ka  =  6.28,  Lossy  layer 


Variable  Frecaiencv  Data 

Figures  16  through  18  are  the  data  for  the  lossless, 
slightly  lossy  and  lossy  -cases  for  b  =  3a  and 
increasingfrequency  for  incident  TM  polarization.  Figure  19 
through  21  are  the  same  data  for  incident  TE  polarization. 
Figure  22  is  the  data  for  the  very  lossy  case  and  incident 
TE  polarization.  The  radius  of  the  conducting  cylinder,  a, 
varies  from  0.003  wavelengths  to  0.4  wavelengths  as  k  is 
increased  from  10  to  1210.  For  the  thickness  of  the  coating 
less  than  approximately  0.4  wavelengths  (kb  <  3.6),  the  two 
solutions  agreed  well  for  both  polarizations.  For  higher 
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Figure  16.  RCS  for  b  =  3a,  Lossless  Layer 


TM  Polarization,  Variable  Frequency 


0  2  4  6  8 

tb 

□  Exact  *  OSfiC 


Figure  17.  RCS  for  b  =  3a,  Slightly  Lossy 

Layer 
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SCS/b  [db]  I  I  ICS/B 


Figure  18.  RCS  for  b  =  3a,  Lossy  Layer 


Figure  19.  RCS  for  b  =  3a,  Lossless  Layer 
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Figure  20.  RCS  for  b  =  3a,  Slightly  Lossy 

Layer 


Figure  21.  RCS  for  b  =  3a,  Lossy  Layer 


56 


Figure  22 


RCS  -for  b  =  3a,  Very  Lossy  Layer 


frequencies,  the  loss  of  the  dielectric  layer  had  to 
increase  for  the  two  solutions  to  agree  over  a  wide  range  of 
frequencies.  For  TM  polarization.  Figure  18,  the  lossy  case 
produced  good  agreement  over  all  the  frequencies  tested.  For 
TE  polarization,  the  very  lossy  case  produced  good  agreement 
over  all  the  frequencies  tested. 

Figures  23  through  25  are  the  data  for  the  lossless, 
slightly  lossy  and  lossy  cases  and  b  =  8.5a  for  incident  TM 
polarization.  For  all  three  cases,  the  OSRC  solution  agreed 


with  the  exact  solution  whenever  kb  was  less  than 
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TH  Polarization,  Variable  Frequency 


□  Exact  ■¥  OSRC 


Figure  25.  RCS  for  b  =  8.5a,  Lossy  Layer 


approximately  3.0  (thickness  of  the  dielectric  layer  less 
than  0.4  wavelengths)  or  kb  became  "very"  large  (a  typical 
high  frequency  approximation) .  For  the  lossless  case, 

Figure  23,  kb  had  to  increase  to  18.5  for  the  OSRC  solution 
to  agree  with  the  exact  solution  again.  For  the  slightly 
lossy  case.  Figure  24,  kb  had  to  increase  to  8.5  for  the 
twosolutions  to  agree.  For  the  lossy  case.  Figure  25,  the 
two  solutions  agreed  for  all  frequencies  tested. 

Figures  26  through  29  are  the  data  for  all  four  cases  and 
b  =  8.5a  for  incident  TE  polarization.  As  in  the  TM 
polarization  data,  for  kb  less  than  3  the  OSRC  solution 
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TE  Polarization,  Variable  Frequency 


ta=0.1,  ta=3.<.  Er=2.54*11.0.  Dr=1.0 


Figure  28.  RCS  for  b  =  8.5a,  Lossy  Layer 


TE  Polarization,  Variable  Frequency 


lia=0.4.  »:b=3.4.  Er=2.S4*11.5,  Di:=1.0 


Figure  29.  RCS  for  b  =  8.5a,  Very  Lossy  Layer 
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produced  good  results.  For  the  lossless  case.  Figure  26,  kb 
less  than  3  was  the  only  region  where  the  OSRC  solution  was 
accurate.  For  the  slightly  lossy  and  lossy  cases.  Figures 
27  and  28,  the  kb  product  had  to  increase  to  9.0  and  6.0, 
respectfully,  for  the  two  solutions  to  agree.  For  the  very 
lossy  case.  Figure  29,  the  two  solutions  agreed  for  all 
frequencies  tested. 

Bistatic  Data 

Figures  30  through  35  are  the  bistatic  data  for  several 
input  conditions.  Figures  30  and  31  are  the  bistatic  RCS 
for  a  "thin",  lossless  dielectric  layer  for  both 
polarizations.  The  input  parameters  were  ka  =  1.57,  kb  = 

3.0  and  lossless  coating.  This  is  equivalent  to  a  coating 
thickness  of  0.23  wavelengths.  There  is  excellent  agreement 
(less  than  1  db  difference)  for  all  angles.  Theta  equal  to 
180  degrees  is  the  monostatic  RCS  point.  Figures  32  and  33 
are  for  a  "thick",  lossless  dielectric  layer  for  both 
polarizations.  The  input  parameters  were  ka  =  1.57  and  kb  = 
6.4.  This  is  equivalent  to  a  coating  thickness  of  0.77 
wavelengths.  The  OSRC  solution  had  significant  error  for 
angles  greater  than  80  degrees  for  TM  polarization  and  100 
degrees  for  TE  polarization.  The  OSRC  solution  was  accurate 
over  a  wide  range  of  angles  in  the  forward  scattering 
direction  for  even  the  thick  coating  case.  Figures  34  and 


TE  Polarization,  Thin  Coating 

113=1.57.  111=3.0,  !r=2.54.  Dr=1.0 
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Figure  31.  Bistatic  RCS,  Thin  Lossless  Layer 
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TM  Polarization.  Thick  Coating 


lca=l.S7.  tli=6,4.  lr=2.S4.  Dr=1.0 


Figure  32.  Bistatic  RCS,  Thick  Lossless  Layer 


TE  Polarization.  Thick  Coating 
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Figure  33.  Bistatic  RCS,  Thick  Lossless  Layer 
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Figure  35.  Bistatic  RCS,  Thick  Lossless  Layer 
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35  are  for  a  "thick"  ,  lossy  layer.  The  input  parameters 
are  the  same  as  for  the  lossless  case  above  except  for  the 
relative  permittivity.  For  TM  polarization,  the  OSRC 
solution  agreed  with  the  exact  solution  for  all  angles.  For 
TE  polarization,  the  OSRC  had  significant  error  for  43°  <  © 

<  130°. 


Data  Analysis 

As  predicted  by  Equation  3.53,  the  behavior  of  the  OSRC 
solution  is  quite  complicated  and  depends  on  many  different 
parameters.  The  magnitudes  of  kb  and  ka,  the  amount  of  loss 
and  thickness  of  the  dielectric  layer  all  played  a  role  in 
determining  the  overall  accuracy  of  the  OSRC  solution.  The 
mechanism  which  produced  this  error  will  be  examined  first, 
followed  by  a  summary  of  the  combinations  of  input 
parameters  which  produced  an  accurate  OSRC  approximation. 

Only  a  few  of  the  terms  of  the  infinite  series  in 
equations  3.49  and  3.50  significantly  contribute  to  the 
magnitude  of  the  overall  sum.  Therefore,  the  error  terms 
presented  in  Equation  3.53  must  remain  small  compared  to  the 
magnitude  of  the  coefficient  for  all  n  less  than  some 
maximum  value  N.  In  order  to  examine  the  behavior  of  the 
error  in  the  OSRC  technique  as  a  function  of  n,  a  computer 
program  was  written  to  compute  the  error  terms  in  Equation 
3.53  for  any  set  of  input  parameters.  Only  first  two  terms 
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of  the  geometric  series  were  used.  A  listing  of  this 
program  is  included  in  Appendix  C.  Figures  36  and  37  show 
the  results  from  the  points  kb  ~  6.4  on  Figure  4  (Figure 
36) ,  a  lossless,  "thick”  dielectric  layer  example,  and 
Figure  6  (Figure  37) ,  a  lossy,  "thick"  dielectric  layer 
example.  These  figures  indicate  that  the  error  is  small  for 
both  the  lossy  and  lossless  cases  for  small  orders,  n  less 
than  6.  However,  for  the  lossless  case,  the  magnitude  of 
the  error  grows  large  compared  to  the  magnitude  of  the 
coefficient  for  larger  orders.  For  example  in  Figure  36, 
the  magnitude  of  the  error  exceeds  tha  magnitude  of  the 
coefficient  for  n  greater  than  7  and  is  almost  half  the 
magnitude  of  the  coeffici-ent  for  n  equal  to  7.  For  the 
lossy  case,  the  magnitude  of  the  error  terms  also  become 
large  compared  to  the  magnitude  of  the  coefficient  for  large 
orders,  but,  the  magnitude  of  the  coefficients  have 
decreased  to  the  point  that  they  do  not  significantly 
contribute  to  the  overall  sum.  Figure  37  illustrates  this 
behavior. 

The  behavior  of  the  magnitude  of  the  error  terms  explain 
why  the  OSRC  technique  is  accurate  for  a  lossy,  "thick" 
dielectric  layer  and  inaccurate  for  a  lossless,  "thick" 
dielectric  layer.  As  the  thickness  increases,  the  number  of 
terms  which  contribute  to.  the  sum,  N,  increases.  Once  N 
gets  to  approximately  7,  the  magnitude  of  the  error  is 
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Coefficient  and  Error  vs  n 


l;a=l.S7,  lcb=fi.4.  Er=2.54,  llr=1.0 


Figure  36.  Error  Mechanism  for  a  Lossless 

Example 


Coefficient  and  Error  vs  n 


Jta=I.S7,  tl)=6.4.  Er=2.54*n  0.  ITr^l.O 


Figure  37.  Error  Mechanism  for  an  Accurate 

Example 
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comparable  to  the  magnitude  of  the  coefficient.  Loss  in  the 
dielectric  layer  improves  this  behavior  because  it  decreases 
the  magnitude  of  the  higher  order  coefficients  and  the 
higher  order  error  terms.  Larger  values  of  kb  also  improved 
the  accuracy  of  the  OSRC  technique. 

Increasing  kb  caused  the  accuracy  to  increase  because  of 
the  e^  term  in  Equation  3.53.  is  proportional  to  (kb) 

where  a  is  at  least  4.  If  Icb  becomes  large  enough,  then  e^ 
dominates  the  behavior  of  the  error  and  the  OSRC  solution  is 
accurate.  Figures  5  and  8  illustrate  this  behavior. 

The  above  discussion  explains  the  mechanism  that  produced 
the  error  in  the  OSRC  technique.  However,  there  were  wide 
ranges  of  input  parameters  where  the  OSRC  technique  produced 
excellent  results.  Figures  38  and  39  outline  the  regions 
where  the  OSRC  technique  was  found  to  be  a  good 
approximation  for  the  RCS  of  a  dielectric  coated  circular 
cylinder  for  both  TM  and  TE  polarizations.  The  minimum  value 
for  kb  is  ka  (a  =  b)  for  any  given  geometry  and  frequency. 
Using  the  data  presented  above  and  extra  data  generated  from 
the  variable  thickness  programs,  the  ranges  of  input 
parameters  which  produced  less  than  2  db  of  error  were 
plotted.  These  regions  were  divided  along  lines  of 
constant  ka  for  varying  kb  and  loss.  Any  combination  of 
geometry,  frequency  and  loss  which  occurs  above  its  line  of 
constant  ka  on  Figure  36  or  37  should  produce  an  OSRC 
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approximation  which  has"  I'ess  than  2  db  of  error.  Based  on 
the  results  presented  in  Figures  36  and  37,  several 
generalizations  about  the  accuracy  of  the  OSRC  technique  are 
possible . 

The  data  presented  above  indicates  that  in  order  for  the 
OSRC  technique  to  agree  with  the  eigenfunction  solution  one 
of  several  possible  conditions  must  be  true.  The  first 
condition  for  a  valid  OSRC  solution  was  a  "thin"  layer.  If 
the  dielectric  layer  was  less  than  about  0.4  wavelengths  for 
incident  TM  polarization,  then  the  OSRC  solution  was  valid. 
If  ka  was  approximately  equal  to  1.57  for  incident  TE 
polarization  and  the  the  -thickness  of  the  dielectric  layer 
was  less  than  0.4  wavelengths,  then  the  OSRC  approximation 
was  good.  Another  condition  for  a  valid  OSRC  solution 
wassignificant  loss  in  the  dielectric  layer.  The  OSRC 
solution  produced  a  good  approximation  to  the  RCS  whenever 
the  imaginary  part  of  the  relative  permittivity  was 
increased  beyond  a  maximum  value.  This  maximum  value  was 
1.0  for  incident  TM  polarization  or  1.5  for  incident  TE 
polarization.  A  third  condition  was  for  a  very  large  kb 
product.  However,  there  was  no  simple  relationship  which 
predicts  how  large  kb  must  be  in  order  for  the  OSRC 
technique  to  have  less  than  2  db  of  error. 


71 


Suminarv 


Comparisons  between  the  OSRC  solution  and  the 
eigenfunction  solution  were  presented.  The  OSRC  solution 
was  found  to  agree  with  the  eigenfunction  solution  for  wide 
ranges  of  input  parameters.  In  general  these  parameters 
were  the  thickness  of  the  dielectric  layer  less  than  0.4 
wavelengths,  the  magnitude  of  the  imaginary  part  of  the 
relative  permittivity  greater  than  1.0  for  incident  TM 
polarization  or  1.5  for  incident  TE  polarization,  or  the  kb 
product  "very"  large.  Data  was  presented  which  outlined 
which  combinations  of  geometry,  dielectric  material  loss  and 
frequency  will  produce  a  valid  OSRC  approximation  for  the 
RCS  of  a  dielectric  coated  cylinder. 


V.  Conclusions  and  Recommendations 


The  OSRC  technique  can  be  used  to  predict  the  RCS  of 
dielectric  coated  cylinders  if  certain  conditions  are  true. 
Figures  38  and  39  summarize  the  regions  of  validity  of  the 
OSRC  technique  for  this  problem.  The  behavior  of  the  OSRC 
solution  was  correct  for  even  surprisingly  small  values  of 
kb.  However,  the  behavior  was  incorrect  for  large  orders. 
This  combination  of  behaviors  made  it  very  difficult  to 
define  a  simple  region  of  validity  based  on  frequency  alone. 
Thus,  Figures  38  and  39  illustrate  which  combinations  of 
loss,  frequency  and  thickness  of  the  dielectric  layer 
produced  accurate  results.  The  data  and  analysis  presented 
in  Chapter  4  are  consistent  with  a  physical  understanding  of 
the  behavior  of  the  field  in  the  dielectric  material. 

Before  summarizing  the  results,  a  brief  discussion  of  these 
physical  mechanisms  will  be  presented. 

The  OSRC  technique  produced  accurate  results  for 
incident  TM  polarization  and  a  thin  layer  because  the 
specular  return  from  the  surface  of  the  cylinder  dominates 
the  RCS  for  this  set  of  conditions.  The  specular  return 
consists  of  the  lowex'  order  modes  (small  indices)  of  the 
infinite  series.  Since  the  OSRC  technique  has  the  correct 
behavior  for  the  lower  order  modes,  there  was  minimal  error 
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in  the  OSRC  solution  for  this  set  of  conditions. 

Incident  TE  polarization  is  much  more  likely  to  excite 
traveling  waves  or  creeping  waves  on  the  surface  of  the 
conducting  cylinder  since'  there  is  a  component  of  the 
electric  field  along  the  surface  in  the  direction  of 
propagation  (Knott  and  others,  1985:166).  These  surface 
waves  are  higher  order  modes  (large  indices)  of  the  infinite 
series  and,  as  explained  above,  the  OSRC  technique  does  not 
behave  well  for  higher  order  modes.  This  explains  why  in 
Figures  13  and  14  the  error  was  large  for  thin  layers; 
creeping  waves  are  contributing  significantly  to  the  RCS  and 
the  OSRC  technique  does  not  approximate  these  higher  order 
modes  well.  Another  mechanism  which  produces  significant 
higher  order  modes  in  the  dielectric  layer  is  multiple 
reflections  between  the-  dielectric  layer  and  the  conducting 
surface . 

Using  geometrical  optics  terminology,  an  incoming  ray  is 
both  transmitted  and  reflected  at  the  dielectric/free  space 
boundary.  The  transmitted  ray  is  then  reflected  at  the 
conductor  surface  and  sent  back  to  the  dielectric  boundary 
where  it  is  both  transmitted  into  free  space  and  reflected 
back  to  the  conductor  surface.  The  reflected  ray  propagates 
back  to  the  conductor  surface  where  it  is  reflected  again 
and  sent  back  to  the  dielectric  boundary.  This  multiple 
reflection  process  continues  with  the  magnitude  of  the 


reflected  ray  decreasing  with  each  reflection  at  the 
dielectric  boundary.  As  the  thickness  of  a  lossless  layer 
increases,  the  contribution  to  the  RCS  from  these  multiple 
reflections  increases.  One  reason  this  occurs  is  simply 
because  the  dielectric  surface  intercepts  a  larger  portion 
of  the  incident  plane  wave.  These  multiple  reflections  are 
higher  order  modes  of  the  infinite  series  and  the  error  in 
the  OSRC  approximation  is  significant  for  these  modes.  If 
there  is  loss  in  the  dielectric  layer,  these  multiple 
reflections  are  attenuated  with  each  pass  through  the  layer, 
thus  reducing  their  contribution  to  the  RCS.  This  physical 
argument  explains  why  the  error  in  the  OSRC  approximation 
gets  large  as  the  thickness  of  the  dielectric  layer 
increases.  These  multiple  reflections  are  higher  order 
modes  of  the  infinite  series  and,  as  explained  many  times 
above,  the  OSRC  approximation  falls  apart  when  the  higher 
order  modes  have  a  significant  contribution  to  the  RCS. 

Loss  in  the  dielectric  layer  reduces  the  contribution  of 
these  higher  order  modes  more  than  the  lower  order  modes  and 
thus  reduces  the  error  in  the  OSRC  approximation.  This 
behavior  is  independent  of  polarization,  and  examination  of 
the  data  in  Chapter  4  is  consistent  with  this  explanation. 
With  the  above  understanding  of  the  physical  mechanisms 
within  the  dielectric  layer  and  the  results  of  the  data 
analysis,  it  is  possible  to  make  some  generalizations  about 
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the  behavior  of  the  OSRC  approximation. 

One  possible  generalization  was  if  the  imaginary  part  of 
the  relative  permittivity  was  greater  the  1.0  for  incident 
TM  polarization  or  1.5  for  incident  TE  polarization,  then 
the  OSRC  technique  produced  a  good  approximation  for  the 
RCS,  even  when  the  size  of  the  cylinder  was  in  the  resonance 
region.  Since  TE  polarization  tends  to  excite  the  higher 
order  modes  in  the  dielectric  layer  more  than  TM 
polarization,  TE  polarization  should  require  more  loss  to 
achieve  the  same  degree  of  accuracy.  This  is  consistent 
with  the  observed  behavio'r.  Many  of  the  materials  used  for 
absorbing  coatings  meet  these  requirements  for  the  imaginary 
part  of  the  relative  permittivity  over  a  wide  range  of 
frequencies  (Joseph,  1986:53-57  and  Eccosorb,  1980:5-20). 

Another  possible  generalization  was  if  the  thickness  of 
the  dielectric  layer  was  less  than  0.4  wavelengths  then  the 
OSRC  approximation  was  accurate  for  TM  polarization.  This 
result  is  very  similar  to  the  range  of  accuracy  found  by 
Wang  for  the  geometrical  optics  approximation  (Wang,  1985: 
962) .  This  similarity  may  be  due  to  the  fact  that  both 
techniques  are  essentially  high  frequency  approximations. 

The  problems  of  numerical  instabilities  near  interior 
resonance  regions  that  occurs  for  MM  solutions  (V'/u  and  Tsai, 
1977:520-523)  were  not  observed.  Whenever  the  input 
parameters  were  such  that  the  error  was  less  than  2  db,  the 


OSRC  solution  was  accurate  even  when  there  was  a  great  deal 
of  cancellation  of  the  fields.  However,  these  deep  null 
regions  were  the  last  areas  to  converge  to  the  exact 
solution  as  the  parameters  changed.  The  OSRC  technique  can 
be  a  valuable  tool  for  cases  where  absorbing  material  has 
been  applied  to  a  conducting  surface. 

The  OSRC  solution  did  not  converge  any  faster  than  the 
eigenfunction  solution.  The  number  of  terms  in  the  infinite 
series  which  had  significant  contributions  to  total  sum  was 
approximately  the  same  for  both  solutions.  Based  on  this 
fact,  the  OSRC  approximation  showed  no  improvement  in 
numerical  efficiency  over  the  exact  solution.  However,  the 
OSRC  technique  is  amenable  to  non-circular  geometries  while 
the  eigenfunction  solution  is  obviously  not. 

Based  on  the  results  presented  above,  there  is 
considerable  motivation  for  continuing  work  with  the  OSRC 
approximation  technique.  The  real  advantage  to  using  this 
technique  should  be  for  non-circular  geometries  where  an 
exact  solution  is  not  known.  Therefore,  the  first 
recommendation  for  future  work  is  to  apply  the  OSRC 
technique  to  an  elliptical  cylinder  coated  with  an 
elliptical  dielectric  layer.  A  numerical  solution  of  the 
interior  problem  using  an  OSRC  boundary  operator  in  lieu  of 
a  radiation  condition  can  be  used  to  find  values  for  the 
fields  on  the  surface  of  the  dielectric  material.  Once 
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these  are  known,  a  quadrature  technique  can  be  used  to  find 
a  numerical  approximation  to  the  RCS  based  on  the  OSRC 
integral  equations  from  Chapter  3. 

Another  recommendation  for  future  work  is  to  extend  the 
present  work  using  OSRC  to  three  dimensions.  A  spheroidal 
shape  would  be  an  excellent  starting  point  for  three 
dimensional  problems  because  the  radius  of  curvature  is 
continuous  over  the  entire  surface.  A  closed  form 
approximation  for  the  RCS  may  be  achievable  using  the  OSRC 
technique. 

A  less  ambitious  project  would  be  to  use  the  programs 
developed  in  this  project  to  investigate  how  well  the  OSRC 
technique  performs  for  magnetic  coatings.  Since  magnetic 
coatings  are  generally  used  to  reduce  traveing  or  creeping 
waves  (Knott  and  others,  1985:253),  the  OSRC  approximation 
may  perform  better  for  this  class  of  coatings. 
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Appendix  A:  Derivation  of  the  OSRC  Operator 
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In  this  Appendix  the  motivation  for  using  an  OSRC 

#  boundary  operator  is  explained  and  the  second  order  operator 
for  a  two  dimensional  field  is  derived. 

Since  the  scattered  field  must  satisfy  the  Sommerfield 

•  radiation  condition,  it  can  be  expressed  as: 


ejkr 


(A.l) 


where, 

k  =  wave  number 

r  =  cylindrical  coordinate  radial  distance 

Equation  A.l  is  valid  for  two  dimensional  fields;  for  three 
dimensional  fields,  the  square  root  of  r  is  replaced  by  r. 
Colton  and  Kress  have  shown  that,  as  a  consequence  of  the 
radiation  condition,  this  series  converges  absolutely  and 
uniformly  and  can  be  differentiated  term  by  term  with 
respect  to  any  of  the  coordinate  system  variables  to  produce 
a  series  which  converges  absolutely  and  uniformly  (Colton 
and  Kress,  1983:72-74).  The  order  OSRC  differential 

operator  uses  this  fact  to  find  an  approximation  foi  the 
normal  derivative  of  the  scattered  field  by  taking  n 


derivatives  of  Equation  A.l  mutliplied  by  e  r^ 

2  n  *“1/2 

error  in  this  approximation  is  proportional  to  r  '  . 


The 
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This  approximation  for  the  normal  derivative  is  used  in 
the  integrand  of  an  integral  equation  representation  for  the 
scattered  field  to  find  an  approximation  for  the  scattered 
field  (Kriegsmann  and  others,  1987:155).  The  details  of  the 
derivation  of  the  second  order  OSRC  operator  follow. 

For  n  =  2,  Equation  A.l  is  multiplied  by  e  j-2/2  . 


-jkr  3/2  ^s 

e  J  X  '  i  = 


^0^  +  ^1  +  ^2  +  ^3  +  .  .  .  . 


Taking  two  derivatives  with  respect  to  r  gives: 


(A. 2) 


e'jkr  r^/2 


0(r) 


(A. 3) 


where  0(r)  ^implies  the  s’eries  behaves  as  r'”  as  r  goes  to 
positive  infinity.  Expanding  the  left  hand  side  of  Equation 


3/2  -jkr  , 
X  '  e  ( 


5^ 


(  3r“  -  2jk)  + 


J  r  2  -  3jkr'^  -  k^  )  )  =  0(r"^  )  (A. 4) 


Since  $  must  satisfy  the  homogenous  Helmholtz  equation,  an 
expression  for  the  second  partial  derivative  of  the 
scattered  field  with  respect  to  r  can  be  found  by  expanding 
the  Laplacian  operator  ia  cylindrical  coordinates  and 
rearranging  terms. 
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6= 


=  -  r 


-1 


S 


-2 


k  ^  (A. 5) 


Sr"  Sr  S  " 

Substituting  Equation  A. 5  into  Equation  A. 4,  dividing  by 
j-3/2  g  3kr  simplifying  the  resulting  expression  produces 


the  second  order  OSRC  operator, 


B2  {  }  =  0(r~'^'^) 


(A. 6) 


where, 


s 

B  =  —  +  ^  -  jk  -  (  — 2  +  ^  )  [  2r(l  -  jkr)]"^ 

Sr  S 


Notice  that  setting  B  {  $  )  =0  will  produce  an  error 

^  s 


term  that  is  proportional  to  r 


•4.5 
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The  algorithms  used  to  compute  the  Bessel's  functions 
were  adapted  from  the  method  described  by  Kelly  (Kelly, 
1967:239-241).  It  is  assumed  that  only  integer  orders  are 
required  and  the  maximum  order  (N)  is  known  a  priori.  The 
algorithms  are  valid  for  either  real  or  complex  argument  and 
a  different  set  of  programs  was  written  for  the  two  cases. 
The  programs  were  written  in  Fortran  77  and  run  on  an  IBM  PC 
machine.  A  listing  of  the  subroutines  is  included.  The 
method  used  to  find  the  Bessel  functions  of  the  first  kind 
(  J^(z)  )  will  be  described  first,  followed  by  a  description 
of  the  method  used  to  find  the  Bessel  functions  of  the 
second  kind  (  ). 

The  Bessel  functions  of  the  first  kind  are  found  using 
backward  racurrsion  and  a  modified  form  of  Gegenbauer's 
formula  to  find  the  normalization  constant.  The  recursion 
relation  is  initialized  by: 


'^N+12 


(z) 


0 


and 


N+11 


(z)  =  10 


r37 


Then  the  initial  values  for  computed  from: 

Jn-i(2>  =  (2n/z)J^(z)  -  (B.l) 

for  1  <  n  <  N+IO.  The  upper  limit  of  N+10  was  chosen  to 
ensure  the  accuracy  of  the  highest  order  Bessel  function 
required  for  use  by  the  calling  program.  These  initial 
values  are  used  in  a  modffied  form  of  Gegenbauer's  formula 
to  find  the  normalization  constant. 
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Kelly  gives  Gegenbauer's  formula  as: 


CO 

2  i2n+a)  r(n+a).  (2/z)“  =  1  (B.2) 

n=0  n! 

where  r(x)  is  the  Gainma  function.  For  n  a  positive  integer 
and  a  =  1,  This  equation  reduces  to: 

n=oo 

(2/z)  E  (2n+l)  J  (z)  =  1  (B.3) 

n=0 

Since  the  satify  the  reccurrence  relation,  they  must 

be  equal  to  a  constant  times  the  Bessel  function. 

J^(z)  =  K  J^(z)  (B.4) 

Substituting  Equation  B.4  into  Equation  B.3,  we  can  solve 
for  the  constant  K. 


K 


(2/z) 


n=oo 

S  (2n+l) 
n=0 


'^2n+i 


(z) 


-1 


(B.5) 


Since  only  a  finite  number  of  Bessel  functions  can  be 
computed,  the  upper  limit  of  the  summation  is  M=(N+10-l)/2 
and  Equation  B.5  is  an  approximation  for  K. 

Once  K  is  found,  each  Jj^(z)  is  multiplied  by  K  for  an 
approximation  to  Jj^(z)  and  Jq(z)  is  found  from  Equation  B.l. 

Once  the  (z)  are  known,  the  (z)  are  found  by  first 
computing  Yq(z)  and  then  using  the  Wronskian  relation  to 
find  the  rest  of  the  Y^(z).  Yq(z)  is  computed  using  the 
following  formula  (Abromawitz  and  Stegun,  1972:360): 
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(z/2)^ 

Yq(z)  -  2/7r[{ln(z/2)+5}J^(z)  +  - 2  "  (1+V2) 


(z/2) 


(1+1/2+1/3) 


(z/2) 

(31) 


ra' 

2 


(1!) 


(2!) 


+ 


] 


(B.6) 


where  5  =  Eulers  constant.  Forward  recurrsion  through  the 
Wronskian  equation  is  used  to  compute  the  rest  of  the 


Y  ( z )  '  s . 
n  ^  ' 


^  "  [^/(Trz)]  }/  Jj^(z)  (B.7) 


The  above  equations  were  implemented  in  the  subroutines 
JNZ,  YNZ  and  YOF  for  complex  argument  and  JNX,  YNX  and  YOFX 
for  real  argument  using  double  precision  arithmetic.  The 
values  for  the  Bessel's  functions  produced  by  these 
subroutines  were  checked  against  the  National  Bureau  of 
Standards  tables  (Abromawitz  and  Stegun,  1972:370-383)  for 

arguments  of  0.5,  2.0  and  10.0  and  orders  0  through  40.  The 

-5 

worst  case  error  found  was  10  when  the  requested  order  was 
at  least  larger  than  the  magnitude  of  the  requested  argument 
plus  three.  A  listing  of  the  programs  follows. 


c 

c 

c 

c 

c 

c 


Subroutine  JNZ ( Z , J , JO , K) 

Computes  integer  order  Bessel  functions  (first  kind)  for 
complex  argument. 

Z  =  argument  of  the  Bessel  Function 
J  =  Complex  valued  array  of  integer  order  Bessel 
functions  K  =  Integer  corotaining  value  of  the  highest 
order  Bessel  function  to  be  computed 


Complex*16  Z , J ( 101) , SUM, F, JO 
Integer  K,LK 


86 


10 


20 


30 


c 

c 

c 


10 


c 

c 

c 

c 

c 

c 

c 

c 


LK=K+10 

IF(float(LK)/2.0.eq.INT(LK/2.0) )  LK=LK+1 
J(LK+2)=(0. 0,0.0) 

J(LK+l)=(1.0E-37,0.0) 
do  10  I  =  LK,1,-1 

J(I)=2.0*Float(I+l)/Z*J(I+l)-J(I+2) 

continue 

sum=J(l)+3.0*J(3) 
do  20  I  =  2, (LK-l)/2 

suin=suin+J  (2. 0*1 +J.)  *(1. 0+2. 0*Float  (I)  ) 
continue 
sum=suiii*  (2 . 0/Z) 
do  30  I  =  1,K 
J(I)=J(I)/suin 
continue 

J0=(2.0/Z)  *J(1)-J(2) 

return 

end 


Subroutine  YNZ (Z, J, J0,Y, Y0,N) 

Complex*16  Z,J(101) ,J0,Y(101) ,Y0,Y0F 

Integer  N 

PI=3. 14159265 

Y0=Y0F(Z, J0,N) 

Y(l)=(Y0*J(l)-2.0/PI/Z)/J0 

Do  10  1=2, N 

Y(I)=(Y(I-l)*J(I)-2.0/PI/Z)/J(I-l) 

Continue 

return 

end 


Coinplex*16  Function  Y0F(Z,J0,N) 

Computes  zeroth  order  Bessel  function  of  the  second 
kind.  Required  input  is  the  complex  valued  argument, 
zeroth  order  Bessel's  function  of  the  first  kind 
evaluated  at  that  argument  and  the  highest  integer  to 
sum  to  (converges  at  1/factorial  squared) 

Complex*16  Z, JO, sum,numer 

Double  precision  K, PI, Eulers, denom, factor 
Integer  N 
PI=3 . 14159265 
Eulers=0. 57721566490153 


c 

c 

c 


initialize  constants  used  in  the  summing  loop 
sum= ( 0 . 0 , 0 . 0 ) 


87 


10 


c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 


10 


20 


30 


c 

c 

c 

c 


numer=Z**2 . 0/4 . 0 
denora=l.  0 
factor=0. 0 
do  10  1=1, N 
K=float(I) 

Factor=Factor+l . 0/K 
denoin=denoin*K 

suin=suin+  (-1 .0)  **(!-!)  *Factor*numer**I/denoiti**2 . 0 
continue 

y0F=2.0/PI*(  (cdlog(Z/2.0)+Eulers)  *J0+suin) 

return 

end 


Subroutine  JNX(X, J, J0,K) 

Computes  integer  order  Bessel  functions  (first  kind) 
for  real  argument. 

X  =  argument  of  the  Bessel  Function 
J  =  Complex  valued  array  of  integer  order  Bessel 
functions  K  =  Integer  containing  value  of  the  highest 
order  Bessel  function  to  be  computed. 


Double  Precision  X, J(lOl) ,SUM, JO 

Integer  K,LK 

LK=K+10 

IF(float(LK)/2.0.eq.INT(LK/2.0) )  LK=LK+1 
J(LK+2)=(0. 0,0.0) 

J(LK+l)=(1.0E-37,0.0) 
do  10  I  =  LK,1,-1 

J(I)=2.0*Float(I+l)/X*J(I+l) -J(I+2) 
continue 

sum=J(l)+3.0*J(3) 
do  20  I  =  2, (LK-l)/2 

sum=sum+J (2.0*1+!)* (1.0+2 . 0* Float (I) ) 
continue 
sum=sum* ( 2 . 0/X) 
do  30  I  =  1,LK 
J(I)=J(I)/sum 
continue  .  . 

J0=(2 . 0/X)  *J(1)  -J(2) 

return 

end 


Double  Precision  Function  Y0FX(X,J0,N) 

Computes  zeroth  order  Bessel  function  of  the  second 
kind. 
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c 

c 

c 

c 

c 


Required  input  is  the  real  valued  argument,  zeroth  order 
Bessel's  function  of  the  first  kind  evaluated  at  that 
argument  and  the  highest  integer  to  sum  to  (converges  at 
1/factorial  squared) 

Double  precision  K, PI , Eulers, denom, numer 
Double  Precision  f actor, X, JO, sum 
Integer  N 
PI=3. 14159265 
Eulers=0. 57721566490153 


c 

c 

c 


10 


initialize  constants  used  in  the  summing  loop 

sum= ( 0 . 0 , 0 . 0 ) 
numer=X**2 . 0/4 . 0 
denom=l. 0 
factor=0 . 0 
do  10  1=1, N 
K=float (I) 

Factor=Factor+l . 0/K 
denom=denom*K 

sum=sum+ (-1 . 0) **(!-!) *Factor*numer**K/denom**2 . 0 
continue 

Y0FX=2 . 0/PI* ( (dlog (X/2 . 0) +Eulers) *J0+sum) 

return 

end 
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Appendix  C.  Listing  of  Computer  Programs 


All  of  the  programs  were  writen  in  FORTRAN  77  and 
implemented  on  an  IBM  PC  computer.  Since  the  code  for  most 
of  the  programs  is  very  simimlar,  only  four  programs  are 
listed  here.  Three  of  the  programs  are  the  programs  used  to 
calculate  the  eigenfunction  solution  and  the  two  OSRC 
solutions  for  the  increasing  thickness  data.  The  fourth 
program  is  the  program  used  to  find  the  error  terms  in  the 
OSRC  technique  as  a  function  of  the  order.  The  only 
difference  between  the  listed  code  and  the  other  programs  is 
the  variable  which  was  changed  after  each  loop. 

The  first  program  is  the  eigenfunction  solution  for  the 
RCS  of  a  dielectric  coated  cylinder  derived  in  Chapter  3. 

C234567890 

Program  DICYLVl 

Complex*16  JNK2B(101) ,J0K2B, JNK2A(101) , J0K2A 
Complex*16  YNK2B(101) , Y0K2B, YNK2A ( 101) ,Y0K2A 
Complex* 16  JNPK2B , YNPK2B, JNPK2A , YNPK2 A , Index , K2 , ANNO 
Complex*16  K2A,K2B,CDSUM, J,DN(2) , AN, ANN, Er ,Ur, Pr 
Double  Precision  PI, JNPM, JN(lOl) , JO, theta 
Double  Precision  YNKIB(IOI) , YOKIB, KIB, YNPKIB 
Double  Precision  Cond,Kl,A,B,Eo,c,DeltaB,X 
Integer  N,L,M,POL 
open ( 8 , f iie= ' outexvl . prn ' ) 
open ( 9 , f ile= ' inv2 . dat ' ) 
c 

c  Initialize  constants 
c 

PI=3. 14159265 
c=3.0E+08  .  . 

Eo=8.854E-12 
J=DCMPLX(0.0, 1.0) 

Read (9,*)  A 

Read (9,*)  B 

Read(9, ' (2D10. 6) ' )  Er 

Read(9, ' (2D10. 6) ' )  Ur 

Read (9,*)  Cond 

Read (9,*)  K1 

Read (9,*)  DeltaB 


Read(9, ' (12) ' )  M 
Read (9,*)  theta 
Read(9, ' (II) ')  POL 
c  POL  =  1  implies  TM  pbrarization 
c  POL  =  0  implies  TE  polarization 
theta=theta*PI/180.0 

Index=CDsqrt (Er*Ur+j*Ur*Cond/ (Kl*Eo*c) ) 

K2=Kl*Index 

K2A=K2*A 

N=int (SORT (Dreal (K2A) **2 . +DIMAG (K2A) **2 . ) +2  0 . ) 
Call  JNZ (K2A, JNK2A, J0K2A,N) 

Call  YNZ(K2A, JNK2A,J0K2A,YNK2A,Y0K2A,N) 

Pr=Ur 

IF(POL.eq.O)  Pr=Er 
c 

c  main  loop 
c 

do  100  L=1,M 
K2B=K2*B 
K1B=K1*B 
B=B+DeltaB 

N=int(SQRT(Dre.al  (K2B)  **2  . +DIMAG  (K2B)  **2. )  )+20. 
c 

c  fill  arrays  with  bessel  functions 
c 

Call  JNX(K1B, JN, J0,-N) 

Call  JNZ(K2B,JNK2B, J0K2B,N) 

Call  YNX(K1B, JN,J0,YNK1B,Y0K1B,N) 

Call  YNZ(K2B,JNK2B, J0K2B,YNK2B,Y0K2B,N) 
If(POL.eq.l)  THEN 

DN(1)=-YNK2B(1) *J0K2A+JNK2B (1) *Y0K2A 
DN (2 ) =Y0K2B*J0K2A-J0K2B*Y0K2A 
ELSE 

DN(1)=YNK2B(1) *JNK2A ( 1) -JNK2B ( 1) *YNK2A(1) 

DN ( 2 ) =- Y0K2  B* JNK2 A ( 1 ) +J0K2  B* YNK2 A ( 1 ) 

ENDIF 

ANN=Pr * JN ( 1 ) *DN ( 2 ) + JO  *Index*DN ( 1 ) 

ANND=ANN+J* (Pr*DN ( 2 ) *YNK1B ( 1) +INDEX*DN ( 1) *Y0K1B) 
CDSUM= J* ANN/ANND* ( - YOKIB* JN ( 1 ) + JO  * YNKl B ( 1 ) ) 
do  40  1=1, N-2 
X=DFLOAT ( I ) 

JNPM=X/K1B*JN(I) -JN(I+1) 

YNPKlB=X/KlB*yNKlB (I) -YNKIB ( I+l ) 
YNPK2B=X/K2B*YNK2B(I)-YNK2B(I+1) 

JNPK2  B=X/K2  B* JNK2  B ( I ) - JNK2  B ( 1+ 1 ) 

IF(POL.EQ.l)  THEN 

DN ( 1 ) = YNPK2  B* JNK2 A ( I ) - JNPK2  B* YNK2 A ( I ) 
DN(2)=YNK2B(I.)  *JNK2A ( I )  -JNK2B  ( I )  *YNK2A(I) 
ELSE 

YNPK2 A=X/K2 A* YNK2 A ( I ) -YNK2 A ( 1+ 1 ) 


JNPK2 A=X/K2 A* JNK2 A ( I ) - JNK2 A ( 1+ 1 ) 

DN ( 1) =YNPK2B*JNPK2A-JNPK2B*YNPK2A 
DN ( 2 ) =YNK2  B ( I ) * JNPK2 A- JNK2  B ( I ) * YNPK2 A 
ENDIF 

ANN=Index* JN ( I ) *DN ( 1) -Pr*JNPM*DN ( 2 ) 

ANND=ANN+J* (Index*YNKlB(I) *DN(1) -Pr*YNPKlB*DN ( 2 ) ) 
AN=J*ANN/MJND*(YNK1B(I)  *JNPM-JN(I)  *YNPK1B) 
CDSUM=CDSUM+2 . 0*AN*dcos (X*theta) 

IF(1.0E05*CDABS(AN) .LT.CDABS(CDSUM) )  goto  41 

40  continue 

41  •  * 

Sigina=PI*KlB**2 . 0/2 . 0*  (real  (CDSUM)  **2 . 0+IMAG  (CDSUM)  **2.0) 
write (8,*)  KIB, Sigma 
100  continue 
stop 
end 


The  next  program  is  the  OSRC  approximation  for  the  RCS  of 
a  dielectric  coated  cylinder  using  the  first  method 
described  in  Chapter  3 . 


C234567890 

Program  OSRCV3 

Complex*16  JNK2B(101) , J0K2B, JNK2A ( 101) ,J0K2A,ANN 
Complex*16  YNK2B(101) ,Y0K2B, YNK2A(101) ,Y0K2A,AN 
Complex*16  JNPK2B,YNPK2B,JNPK2A,YNPK2A, Index, K2 , ANNO 
Complex*16  K2A, K2B, CDSUM, J, DN (3 ) , KlBl , DN30 , Er , Ur , Pr 
Double  Precision  PI, J0,X,theta, JNC(lOl) , JOC 
Double  Precision  JlfPM, JN(lOl)  , YNKIB(IOI)  , Y0K1B,K1B 
Double  Precision  Cond,Kl, A, B, Eo, c, DeltaB, YOFX 
Integer  N,L,M,POL 
open ( 8 , f ile= ' outv3 . prn ' ) 
open ( 9 , f ile= ' inv2 . dat • ) 
c 

c  Initialize  constants 
c 

PI=3. 14159265 
c=3 . OE+08 
Eo=8.854E-12 
J=DCMPLX(0. 0,1.0) 

Read (9,*)  A 

Read (9,*)  B 

Read(9, ' (2D10.6) ' )  Er 

Read (9, ' (2D10.6) ' )  Ur 

Read (9,*)  Cond 

Read (9,*)  K1 


c 

c 


c 

c 


c 

c 

c 


Read (9,*)  Deltas 
Read(9, ' (12) ' )  M 
Read ( 9 , * )  theta 
Read(9, ' (II) ')  POL 
POL  =  1  implies  TM  polarization 
POL  =  0  implies  TE  polarization 
theta=theta*PI/180 . 0 

Index=CDsqrt (Er*Ur+J*Ur*Cond/ (Kl*Eo*c) ) 

K2=Kl*Index 

K2A=K2*A 

N=int(SQRT(real(K2A)**2.+IMAG(K2A) **2. ) )+10. 
Call  JNZ (K2A, JNK2A, J0K2A,N) 

Call  YNZ (K2A, JNK2A, J0K2A, YNK2A, Y0K2A,N) 

Pr=Ur 

IF(POL.EQ.O)  Pr=Er 

main  loop  '  * 

do  100  L=1,M 
K2B=K2*B 
K1B=K1*B 
B=B+DeltaB 
K1B1=1. 0-J*KlB 

N=int(SQRT(real(K2B) **2.+IMAG(K2B) **2. ) )+20. 

fill  arrays  with  bessel  functions 

Call  JNX(K1B, JN, J0,N) 

Call  JNZ (K2B, JNK2B, J0K2B,N) 

Call  YNX(K1B, JN, J0,YNK1B, Y0K1B,N) 

Call  YNZ (K2B, JNK2B, J0K2B, YNK2B, Y0K2B,N) 

IF  (POL.EQ.l)  THEN 

DN(1)=-YNK2B(1) *J0K2A+JNK2B ( 1) *Y0K2A 
DN ( 2 ) =Y0K2B* J0K2A-J0K2B*Y0K2A 
ELSE 

DN(1)=YNK2B(1) *JNK2A(1) -JNK2B(1) *YNK2A(1) 

DN ( 2 ) =-Y0K2B*JNK2A (1) +J0K2B*YNK2A (1) 

ENDIF 

DN30=J-1 . 0/ (2 . 0*K1B)  +1. 0/ (8 . 0*K1B*K1B1) 
ANN=(-JN(1) -DN30*J0) **2.*DN(2) 
ANND=INDEX*DN(1) -Pr*DN(2) *DN30 
CDSUM=ANN/ANND 
do  40  1=1, N-2 
X=Df.oat  (I) 

JNPM=X/K1B*JN(I) -JN(I+1) 

YNPK2B=X/K2B* YNK2B ( I ) -YNK2B ( I+l ) 
JNPK2B=X/K2B*JNK2B(I)-JNK2B(I+1) 
IF(POL.EQ.l)  THEN 

DN ( 1 ) =YNPK2  B* JNK2A ( I ) - JNPK2  B* YNK2 A ( I ) 
DN(2)=YNK2B(I) *JNK2A ( I ) -JNK2B ( I ) * YNK2A ( I ) 
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ELSE 

YNPK2 A=X/K2 A* YNK2 A ( I ) -YNK2 A ( I+l ) 

JNPK2 A=X/K2 A* JNK2 A ( I ) - JNK2 A ( 1+ 1 ) 

DN ( 1) =YNPK2B*JNPK2A-JNPK2B*YNPK2A 
DN(2)=YNK2B(I) * JNPK2A-JNK2B ( I ) *YNPK2A 
ENDIF 

DN(3)=DN30-(X**2 .0)/(2.0*KlB*KlBl) 

ANN=(JNPM-DN(3) *JN(I) ) **2.*DN(2) 

ANND=INDEX*DN ( 1 ) -Pr *DN ( 2 ) *DN ( 3 ) 

AN=ANN/ANND 

CDSUM=CDSUM+2 . 0*AN*Dcos (X*theta) 

IF(1.0E05*CABS(AN) .LT.CABS(CDSUM) )  goto  41 

40  continue 

41  X=2 . 0*dsin( theta/2 . 0) *K1B 
M=INT(X)+3. 

Call  JNX(X, JNC, J0C,M) 

CDSUM=Pr*CDSUM+ (DN30-K1B/  (4 . 0*K1B1) ) *J0C+ (1 . 0+1/ (2 . 0*K1B1) ) 
x*sin ( theta/2. ) *JNC (1) +K1B/ (4 . 0*K1B1) *cos (theta) *JNC(2) 

Sigina=PI*KlB**2 . 0/2 . 0*  (Dreal  (CDSUM)  **2 . 0+DIMAG  (CDSUM)  **2 .0) 
write(8,*)  KIB, Sigma 
100  continue 
stop 
end 


The  next  program  is  the  OSRC  approximation  for  the 
RCS  of  a  dielectric  coated  cylinder  using  the  second  method 
described  in  Chapter  3 . 


C234567890 

Program  OSRCV5 

Complex*16  JNK2B(101) , J0K2B, JNK2A(101) ,J0K2A,AN 
Complex*16  YNK2B(lcri)  ,  Y0K2B,  YNK2A ( 101)  ,Y0K2A,ANND 
Complex*16  JNPK2B,YNPK2B, JNPK2A,YNPK2A, Index, K2 
Complex*16  K2A,K2B, CDSUM, J,DN(3) ,KlBl,Er,Ur, Pr 
Complex*16  ANN,DN30 
Double  Precision  PI, JO, X, theta 

Double  Precision  JNPM, JN (101) , YNKIB ( 101) , YOKIB, KIB 

Double  Precision  Cond, K1 , A, B, Eo, c, DeltaB 

Integer  N,L,M,POL 

open (8, file='outv5.prn' ) 

open (9 , f ile= ' inv2 .dat ' ) 


c  Initialize  constants 
c 

PI=3 . 14159265 
c=3 . OE+08 
Eo=8.854E-12 
J=DCMPLX(0.0, 1.0) 

Read (9,*)  A 
Read (9,*)  B 
Read(9, ' (2D10.6) ' )  Er 
Read(9, • (2D10.6) ' )  Ur 
Read (9,*)  Cond 
Read (9,*)  K1 
Read (9,*)  DeltaB 
Read(9, ' (12) ' )  M 
Read (9,*)  theta 
Read(9, ' (II) ' )  POL 
c  POL  =  1  implies  TM  polarization 
c  POL  =  0  implies  TE  polarization 
theta=theta*PI/180 . 0 

Index=CDsqrt  (Er*Ur+-J*Ur*Cond/  (Kl*Eo*c)  ) 

K2=Kl*Index 

K2A=K2*A 

N=int(SQRT(real (K2A) **2.+IMAG(K2A) **2. ) )+20 
Call  JNZ (K2A, JNK2A, J0K2A,N) 

Call  YNZ (K2A, JNK2A, J0K2A, YNK2A, Y0K2A,N) 
Pr=Ur 

IF(PCL.EQ.O)  Pr=Er 
C 

c  main  loop 
c 

do  100  L=1,M 

K2B=K2*B 

K1B=K1*B 

B=B+DeltaB 

K1B1=1.0-J*K1B 

N^int (SQRT(real (K2B) **2 . +IMAG (K2B) **2. ) )+20 
c 

c  fill  arrays  with  bessel  functions 
c 

Call  JNX(K1B, JN, J0,N) 

Call  JNZ (K2B, JNK2B, J0K2B,N) 

Call  YNX(K1B,  JN,  J0,.YNK1B,  Y0K1B,N) 

Call  YNZ (K2B, JNK2B, J0K2B, YNK2B, Y0K2B,N) 

IF  (POL.EQ.l)  THEN 

DN(1)=-YNK2B(1) *J0K2A+JNK2B ( 1 ) *Y0K2A 
DN(2)=YOK2B*JOK2A-JOK2B*YOK2A 
ELSE 

DN(1)=YNK2B(1) *JNK2A ( 1) -JNK2B ( 1 ) *YNK2A(1) 
DN ( 2 ) =-Y0K2B*JNK2A (1) +J0K2B*YNK2A ( 1 ) 

ENDIF 


n 


DN30=J-1 .0/(2. 0*K1B) +1 . 0/ (8 . 0*K1B*K1B1) 

ANN=(-JN(1) -DN30*J0) *(INDEX*DN(1) * J0+Pr*DN ( 2 ) *JN ( 1 ) ) 
ANND=INDEX*DN(1) -Pr*DN(2) *DN30 
CDSUM=ANN/ANND 
do  40  1=1, N-2 
X=DFloat (I) 

JNPM=X/K1B*JN(I) -JN(I+1) 

YNPK2  B=X/K2  B* YNK2  B ( I ) - YNK2  B ( 1+ 1 ) 

JNPK2  B=X/K2  B* JNK2  B ( I ) - JNK2  B ( 1+ 1 ) 

IF(POL.EQ.l)  THEN 

DN ( 1 ) =YNPK2B* JNK2A ( I ) -JNPK2B* YNK2 A ( I ) 
DN(2)=YNK2B(I) * JNK2A( I ) -JNK2B ( I ) *YNK2A ( I ) 

ELSE 

YNPK2 A=X/K2 A*YNK2A ( I ) -YNK2A ( I+l ) 
JNPK2A=X/K2A*JNK2A(I)-JNK2A(I+1) 

DN ( 1) =YNPK2B*JNPK2A-JNPK2B*YNPK2A 
DN(2)=YNK2B(I) *JNPK2A-JNK2B(I) *YNPK2A 
ENDIF 

DN(3) =DN30- (X**2 . 0)/ (2 . 0*K1B*K1B1) 

ANN={JNPM-DN(3) *JN(I) ) * ( INDEX*DN ( 1) *JN ( I ) -Pr*DN ( 2 ) *JNPM) 
ANND=INDEX*DN(1) -Pr*DN(2) *DN(3) 

AN=ANN/ANND 

CDSUM=CDSUM+2 . 0*AN*cos (X*theta) 

IF(1.0F05*CABS(AN) .LT.CABS(CDSUM) )  goto  41 

40  continue 

41  continue 

Sigina=PI*KlB**2 . 0/2 . 0*  (real  (CDSUM)  **2 . 0+IMAG  (CDSUM)  **2.0) 
write(8,*)  KIB, Sigma 

100  continue 
stop 
end 


The  next  program  was  used  to  calculate  the  error  terms 

and  coefficients  of  Equation  3.53  and  Figures  36  and  37. 

C234567890 

Program  ERROR 

Complex*16  JNK2B(101) , J0K2B, JNK2A ( 101) ,J0K2A,BN(3) , YFC 
Complex*16  YNK2B(101) , Y0K2B, YNK2A ( 101) ,Y0K2A,XN 
Complex*16  JNPK2B, YNPK2B,JNPK2A, YNPK2A, Index, K2 , ANND 
Complex*16  K2A,K2B, CDSUM, J,DN(3) , KlBl , DN3 0 , Er , Ur , Pr 
Complex*  16  YN ,  YNN-,  Y-NND ,  AN ,  AN ,  ANN 

Double  Precision  PI , JO , X, theta , JNC ( 101) , JOC , YOFX 
Double  Precision  JNPM, JN(lOl) , YNKIB(IOI) , Y0K1B,K1B 
Double  Precision  Cond,Kl,A,B,Eo,c,DeltaB 


c 

c 

c 


c 

c 


Integer  N,L,M,POL 
open (8 , file- ' outv3 .prn' ) 
open (9, file=' inv2 .dat' ) 
open (10, f ile= 'EN.prn ' ) 
open (11, file='ENYN.prn' ) 
open  ( 12  ,  f  ile= '  ENXN  ."pm ' ) 
open (14, file=' AN.prn' ) 
open (15, file='XNYN.prn' ) 

Initialize  constants 

PI=3. 14159265 
c=3 . OE+08 
Eo=8.854E-12 
J=DCMPLX(0. 0, 1. 0) 

Read (9,*)  A 
Read(9,*)  B 
Read(9, ' (2D10. 6) ' )  Er 
Read(9, ' (2D10.6) ')  Ur 
Read (9,*)  Cond 
Read (9,*)  K1 
Read (9,*)  DeltaB 
Read(9, ' (12) ' )  M 
Read (9,*)  theta 
Read(9, ' (II) ' )  POL 
POL  =  1  implies  TM  polarization 
POL  =  0  implies  TE  polarization 
theta=theta*PI/180  .-0 

Index=CDsqrt (Er*Ur+J*Ur*Cond/ (Kl*Eo*c) ) 

K2=Kl*Index 

K2A=K2*A 

N=int (SQRT(Dreal (K2A) **2 . +DIMAG (K2A) **2.)+20. ) 
Call  JNZ (K2A, JNK2A,J0K2A,N) 

Call  YNZ (K2A, JNK2A, J0K2A,YNK2A, Y0K2A,N) 

Pr=Ur 

IF(POL.eq.O)  Pr-Er 


c 

c  main  loop 
c 


do  100  L=1,M 
K2B=K2*B 
K1B=K1*B 
K1B1=1. 0-J*KlB 


B=B+DeltaB 

N=int(SQRT(Dreal (K2B) **2,+DIMAG(K2B) **2. ) )+20. 
c 

c  fill  arrays  with  bessel  functions 
c 


Call  JNX(K1B, JN, J0,N) 

Call  JNZ (K2B, JNK2B,J0K2B,N) 


98 


Call  YNX(K1B, JN, JO, YNK1B,Y0K1B,N) 

Call  YNZ(K2B, JNK?B, J0K2B,yNK2B, Y0K2B,N) 
If(POL.eq.l)  TH2N 

DN(1)=-YNK2B(1) *J0K2A+JNK2B ( 1) *Y0K2A 
DN ( 2 ) =Y0K2B*J0K2A-J0K2B*Y0K2A 
ELSE 

DN(1)=YNK2B(1) * JNK2A ( 1) -JNK2B ( 1) *YNK2A(1) 

DN ( 2 ) =-Y0K2B* JNK2 A (1) +J0K2B* YNK2A ( 1 ) 

ENDIF 

DN30=J-1. 0/ (2 . 0*K1B) +1. 0/ (8 . 0*K1B*K1B1) 
BN(1)=(-JN(1)-J*YNK1B(1) ) 

BN(2)=BN(1)/(J0+J*Y0K1B) 

BN(3)=BN(2) -DN30 

Xn=(BN(3) *Pr*DN(2) ) / ( INDEX*DN ( 1) -Pr*DN ( 2 ) *BN ( 2 ) ) 
ANN=Pr*JN (1) *DN ( 2 ) +JO*Index*DN (1) 

ANND=ANN+J* (Pr*DN ( 2 ) *YNK1B ( 1) +INDEX*DN ( 1 ) *Y0K1B) 
CDSUM=J*ANN/ANND* (-Y0K1B*JN ( 1) +J0*YNK1B ( 1 ) ) 
YNN=-A1JN*BN(2)  *JN(1)  -Pr*DN ( 2 )  *CDSUM 
YND=ANN* {J0+BN(2) *JN(1) ) 

YN=YNN/YND  '  * 

write (12,*)  XN*BN(3) *CDSUM 
write (10,*)  BN(3)*CDSUM 
write (11,*)  YN*BN(3) *CDSUM 
write (14,*)  cdsum 
write (15,*)  XN*YN*CDSUM*BN(3) 
do  40  1=1, N-2 
X=DFLOAT(I) 

JNPM=X/K1B*JN(I)-JN(I+1) 

YNPK1B=X/K1B*YNK1B(I) -YNK1B(I+1) 

YNPK2  B=X/K2  B* YNK2  B ( I ) -YNK2  B ( 1+ 1 ) 

JNPK2  B=X/K2  B* JNK2  B ( I ) - JNK2  B ( 1+ 1 ) 

IF(POL.EQ.l)  THEN 

DN ( 1 ) =YNPK2  B* JNK2 A ( I ) -JNPK2  B* YNK2 A ( I ) 
DN(2)=YNK2B(I) *JNK2A(I) -JNK2B(I) *YNK2A(I) 

ELSE 

YNPK2 A=X/K2 A* YNK2 A ( I ) -YNK2 A ( I+l ) 

JNPK2 A=X/K2 A* JNK2 A ( I ) -JNK2 A ( I+l ) 

DN ( 1) =YNPK2B*JNPK2A-JNPK2B*YNPK2A 
DN(2)=YNK2B(I) *JNPK2A-JNK2B(I) *YNPK2A 
ENDIF 

DN(3)=DN30-(X**2.0)/(2.0*K1B*K1B1) 

BN(1)  =  (JNPM+J*YNPK1-B) 
BN(2)=BN(l)/(Jn(I)+J*YNKlB(I) ) 

BN(3)=BN(2) -DN(3) 

Xn=(BN(3) *Pr*DN(2) )/ (INDEX*DN(1) -Pr*DN ( 2 ) *BN (2 ) ) 
ANN=Index*JN(I) *DN(1) -Pr*JNPM*DN (2) 

AN1!D=ANN+J*  (Index*YNKlB (I)  *DN(1)  -Pr*YNPKlB*DN(2)  ) 
AN=o*ANN/ANND* ( YNKIB ( I ) *JNPM-JN ( I ) *YNPK1B) 

YNN= -ANN*BN ( 2 ) * JN ( 1 ) -Pr*DN ( 2 ) *AN 
YND=ANN* ( JN ( I ) -BN ( 2 ) * JNPM) 


99 


yN=YNN/yND 

write (12,*)  XN*BN(3)*AN 
write(10,*)  BN(3)*AN 
write  (11,*)  yN*BN(3J*AN 
write (14,*)  AN 
write (15,*)  XN*yN*AN*BN(3) 

CDSUM=CDSUM+2 . 0*AN*dcos (X*theta) 

IF(1.0E05*CDABS(AN) .LT.CDABS(CDSUM) )  goto  41 

40  continue 

41 

Sigina=PI*KlB**2 . 0/2 . 0*  (real  (CDSUM)  **2 . 0+IMAG  (CDSUM)  **2.0) 
write (8,*)  KIB, Sigma 

100  continue 
stop 
end 


Bibliography 


Abromowitz,  Milton  and  Stegun,  Irene  A.  Handbook  of 
Mathematical  Functions.  Dover  Publications  Inc.,  New 
York.  1965. 

Baker,  William.  Class  Notes  from  MATH  5.12,  Mathematical 
Methods  of  Scattering.  School  of  Engineering,  Air  Force 
Institute  of  Technology  (AU) ,  Wright-Patterson  AFB  OH, 
January  1988. 

Bowman,  J.J.  et  al .  .  Eds.*  Electromagnetic  and  Acoustic 
Scattering  by  Simple  Shapes.  Amsterdam:  North-Holland 
Publishing  Company.  1969. 

Chang,  Shu-Kong  and  Kenneth  K.  Mei.  "Application  of  the 
Unimoment  Method  to  Electromagnetic  Scattering  of 
Dielectric  Cylinders,"  IEEE  Transactions  on  Antennas  and 
Propagation.  25:  35-42  (January  1976) 

Colton,  David  and  Kress,  Rainer.  Integral  Eguations  Methods 
in  Scattering  Theory.  Wiley  and  Sons,  New  York.  1983. 

Gradshteyn,  I.  S.  and  Ryzhik,  I.  M.  Table  of  Integrals, 

Series,  and  Products.  Academic  Press,  New  York.  1980. 

Harrington  Roger  F.  Time-Harmonic  Electromagnetic  Fields. 
McGraw  Hill  Book  Co.,  New  York.  1961. 

-  High  Loss  Dielectrics  Microwave  Absorbers.  Eccosorb 

Product  Catalog,  Emerson  and  Cumriings,  Canton  MA.  1980. 

James,  Graeme  L.  Geometrical  Theory  of  Diffraction  for 
Electromagnetic  Waves.'  Stevenage,  England,  Peter 
Peregrinus  Ltd.  1976. 

Jin,  Jian-Ming  and  Valdis  V.  Liepa.  "Application  of  Hybrid 
Finite  Element  Method  to  Electromagnetic  Scattering  from 
Coated  Cylinders,"  IEEE  Transactions  on  Antennas  and 
Propagation,  36;  50-54  (January  1988) 

Joseph,  Lt  John.  Multiple  Angles  of  Incidence  Measurement 
Technigue  for  the  Permittivity  and  Permeability  of  Lossy 
Material  at  Millimeter  Wavelengths.  MS  Thesis 
AFIT/GE/ENG/86D-58 .  School  of  Engineering,  Air  Force 
Institute  of  Technology  (AU) ,  Wright  Patterson  AFB  OH, 
December  1986. 


Jost  Randy  J.  Class  Notes  from  EENG  6.30,  Applications  of 
Electromagnetic  Theory.  School  of  Engineering,  Air  Force 
Institute  of  Technology  (AU) ,  Wright-Patterson  AFB  OH, 
January  1988. 

Kelly,  Louis  G.  Handbook  of  Numerical  Methods  and 

Applications .  Addison-Wesley  Co.  Reading,  MA.  1967. 

Knott  Eugene  F.  and  others.  Radar  Cross  Section.  Artech 
House,  Inc.  Norwood  MA.  1985. 

Kriegsman,  Gregory  A.  and  Cathleen  S.  Morawetz.  "Solving 
the  Helmholtz  Equation  for  Exterior  Problems  with 
Variable  Index  of  Refraction  :  I,"  SIAM  Journal  for 
Scientific  and  Statistical  Computation.  1:  371-385 
(September  1980) 

Kriegsman,  Gregory  A.  et  al .  "A  New  Formulation  of 

Electromagnetic  Scattering  Using  an  On-Surface  Radiation 
Boundary  Condition  Approach,"  IEEE  Transactions  on 
Antennas  and  Propagation.  35;  153-161  (February  1987) 

Leader  Carl  J.  "Modern  Methods  of  Scattering  Predictions," 
Proceedings  of  SPIE  -  The  International  Society  for 
Optical  Engineering.  Applications  of  Mathematics  in 
Modern  Optics.  358:  17-23  (August  1982) 

Peters,  Leon  Jr.  "End-Fire  Echo  Area  of  Long  Thin  Bodies," 
IRE  Transactions  on  Antennas  and  Propagation.  1;  133-139 
(January  1958) 

Ruck  G.  T.  et  al .  Radar  Cross  Section  Handbook.  Plenum 
Press,  New  York.  1970. 

Tang,  Charles  C.  H.  "Backscattering  from  Dielectric-Coated 
Infinite  Cylindrical  Obstacles,"  Journal  of  Applied 
Physics,  28:  628-633  (May  1957) 

Wang,  Nan.  "Electromagnetic  Scattering  from  a  Dielectric- 
Coated  Circular  Cylinder,"  IEEE  Transactions  on  Antennas 
and  Propagation,  33:  '960-963  (September  1985) 

Wu,  Te-Hao  and  Leonard  L.Tsai.  "Scattering  by  Arbitrarily 
Cross-Sectioned  Layered  Dielectric  Cylinders,"  IEEE 
Transactions  on  Antennas  and  Propagation,  25;  518-524 
(July  1977) 


VITA 


Captain  William  Getting  wa?  born  <^i'^l||imHllimil||||PE^ 

He  graduated  from  high  school  in 
1974.  the 

Force  in  1978  and  was  selected  to  attended  the  University  of 
Washington,  Seattle,  Washington  under  the  Airmans  Education 
and  Commissioning  Program.  He  graduated  Magna  Cum  Laude  in 
December  1983  with  a  Bachelor  of  Science  in  Electrical 
Engineering  and  was  commissioned  in  April  1984  after 
completing  Officers  Training  School  at  Lackland  AFB,  Texas. 
Before  entering  the  School  of  Engineering,  Air  force 
Institute  of  Technology,  he  worked  for  three  years  at  the 
Air  Force  Wright  Aeronautical  Laboratories,  Avionics 
Laboratory  as  a  electronic  warfare  research  engineer. 


UNCLiiooIFIiiu 


SECURITY  CLASSIFICATION  OF  THIS  PAGE 


REPORT  DOCUMENTATION  PAGE 


lb.  RESTRICTIVE  MARKINGS 


Form  Approved 
0MB  No.  0704-0188 


la.  REPORT^SECURITY  CLASSIFICATION 

U  LN  C  i.  A  O  O  X  r’ i  jAU 


2a.  SECURITY  CLASSIFICATION  AUTHORITY 


2b.  DECLASSIFICATION /DOWNGRADING  SCHEDULE 


4.  PERFORMING  ORGANIZATION  REPORT  NUMBER(S) 

AFIi'/ivlA/G£/88D-2 


Y'**.  NAME  OF  PERFORMING  ORGANIZATION  6b  OFFICE  SYMBOL 

i  (lf  applicable) 

xiChool  of  engine uring  AFIf/jAliA 


6c.  ADDRESS  (Gty,  State,  and  ZIP  Code) 

Air  Force  institute  of  'technology 
V/right-Patterson  Ohio  45^33 


8b.  OFFICE  SYMBOL 
(If  applicable) 

wF'V/aL/ Cuj 


3.  DISTRIBUTION /AVAILABILITY  OF  REPORT 

Approved  for  Public  Release; 
distribution  Unlimited 


S.  MONITORING  ORGANIZATION  REPORT  NUMBER(S) 


7a.  NAME  OF  MONITORING  ORGANIZATION 


8c.  ADDRESS  (C/fy,  State,  and  ZIP  Code) 


9.  PROCUREMENT  INSTRUHE 


10.  SOURCE  OF  FUNDING  NUMBERS 


^  . 

PROGRAM 

PROJECT 

TASK 

WORK  UNIT 

•■hBaFo  ,  Ohio  45^33 

ELEMENT  NO. 

NO. 

NO 

ACCESSION  NO. 

11.  TITLE  (Include  Security  Classification) 

ihe  Radar  Cross  Section  of  a  two  dimensional  uielectric  Coated  Cylinder 
Using  0n-8urface  Radiation  Conditions _ UNCLASSIFIRd _ 


12.  PERSONAL  AUTHOR{S) 

William  D.  Getting,  Capt.  USaF  _ 


13a.  TYPE  OF  REPORT  13b.  TIME  COVERED  14.  DATE  OF  REPORT  {Year,  Month,  Day) 

WR  Thesis  from _ to _  I988  December 


16.  SUPPLEMENTARY  NOTATION 


ns 


COSATI  CODES 


GROUP  SUB-GROUP 


18.  SUBJECT  TERMS  {Continue  on  reverse  if  necessary  and  identify  by  block  number) 

Radar  Cross  oections,  xilectromagnetic  Scattering 
islectromagnetic  Theory,  Scattering  Cross 
Sections 


19.  ABSTRACT  {Continue  on  reverse  if  necessary  and  identify  by  block  number) 

Thesis  Chairmans  William  Baker,  Lt  Col,  UbaF 

Assistcint  Professor  of  Applied  Mathematics 

Abstract  s 

The  RCo  of  a  dielectric  coated  circular  cylinder  was  found  using  an  On- 
L/urface  Radiation  Condition  (OSHC)  and  compared  tea  eigenfunction  solution. 

The  OonC  approximation  was  found  to  agree  with  the  eigenfunction  solution 
whenever  the  higher  order,  modes  of  the  field  in  the  dielectric  layer  were 
insignificant.  An  electrically  thin  layer  or  a  resonable  amount  of  loss  in 
the  dielectric  material  would  produce  an  accurate  OoRC  approximation  for 
resonant  size  scatterers.  Specific  results  indidated  that  the  magnitude 

"the  relative  permittivity  should  be  greater  than 
i.U  for  an  accurate  approximation. 


20.  DISTRIBUTION /AVAILABILITY  ABSTRACT 
□  UNCLASSIFIED/UNLIMITED  ©  SAME  AS  RPT. 


Si 


,  21.  ABSTRACT  SECURITY  CLASSIFICATION 

□  DTIC  USERS  UNuLARSIFllLiJ 


22b.TELEPHpNE//nc/ude  Area  Code)  22c.  OFFICE  SYMBOL 


DD  Form  1473,  JUN  86 


22b.J^L^PHgE^,n^.^dgA 


Previous  editions  are  obsolete.  _ SE 


AFIT/ima 


SECURITY  CLASSIFICATION  OF  THIS  PAGE 

UN  C  LAS  R I  FI  RD 


